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Steam in the Ring Discharge 
By G. I. FINCH* 
Presidential Address, delivered 6th May 1949 


ABSTRACT. The behaviour of steam and its decomposition products in the ring 
discharge has been examined. 

Dry hydrogen is not dissociated. "The production of atomic hydrogen is dependent 
upon the presence of steam which dissociates into hydroxyl and atomic hydrogen. A second- 
ary source of atomic hydrogen is then afforded by the interaction of hydroxyl with molecular 
hydrogen. 

The escape from the discharge of atomic hydrogen, a long-lived species, favours the 
dissociation of steam. Mercury vapour, on the other hand, inhibits the formation of atomic 
hydrogen and thus leads to a high equilibrium steam concentration. 

Unlike dry hydrogen, dry oxygen is dissociated into atoms, but these have a short life as 
such and recombine in the discharge to form molecular oxygen and ozone. 

The reaction mechanisms occurring in the discharge are discussed in the light of 
spectrographic results. . 


§1. INTRODUCTION 


ETWEEN 1926 and 1935 my collaborators and I carried out a systematic 

study of the combustion of hydrogen and carbon monoxide by oxygen in 

the cathode glow of the electric discharge maintained between metallic 
electrodes. Our object was to elucidate the nature of the mechanisms of these 
combustion processes, and the results and conclusions drawn therefrom were 
summarized in the Fournal of the Chemical Society (Finch 1935). It may be 
remarked that, during the progress of this work, it became apparent that the re-. 
sults were also of practical interest, in that they threw fresh light on the chemistry 
of the thermionic valve, particularly in relation to the electron emissivity of the 
cathode. 

Whilst certain phenomena arising out of the use of metallic electrodes in our: 
experiments had helped towards an understanding of some aspects of the general 
mechanism of combustion, it became evident that, in other directions, they tended 
to obscure the issues. Indeed, as regards the undoubtedly complex mechanism 
‘of the homogeneous combustion of hydrogen by oxygen to steam, we could, in 
1935, say no more than that, when the combustible mixture was rigidly dried, 
the initial step was the formation of the hydroxyl radical and that, in the undried 
mixture, steam was first dissociated into atomic hydrogen and hydroxyl, whereupon. 
the hydroxyl somehow acted as a powerful promoter in the succeeding stages, 
leading ultimately to the formation of steam. No positive réle in the combustion 

* The Presidential Address will also appear in Proc. Phys. Soc. for September 1949, Section B. 
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process could be assigned to atomic hydrogen. ‘This was remarkable since, under 
the conditions of our experiments, atomic hydrogen was known to be produced 
in considerable quantities. Accordingly it was decided to examine the equili- 
brium between steam and its dissociation products in the electrodeless discharge 
under conditions which would eliminate surface effects, or at least permit of an 
adequate appreciation of the part played by them. 


§2. APPARATUS AND GENERAL PROCEDURE 


The gas (steam, 2H, +O,, Hg, O, or rare gas, or various mixtures thereof) 
was contained in the closed circulation system shown in Figure 1. This was 
assembled from the cleaned components using a dry, filtered, glass-blowing air 
supply, and was so constructed that the circulating gas came into contact only 
with glass, mercury or quartz, thus avoiding the contaminations introduced by 


1) 
QUARTZ | 
— 


RARE G, 


TO HYVAC = RESERVOIR 
HYDROGEN 
OR OXYGEN 


MANOMETER 


ROUGH 
PUMP. 
TO AIR<= 


Figure 1. 


lubricated stopcocks and rubber joints which have in the past so often vitiated the 
results of work on the combustion of hydrogen. ‘The glass and quartz components 
were cleaned with chromic acid and distilled water ; the mercury, previously 
washed in the usual manner and distilled, was redistilled into the appropriate 
vessels ; and the phosphoric oxide was redistilled in oxygen until free from the 
lower oxides (Finch and Fraser 1926). ‘The major remaining sources of impurity 
(other than mercury vapour) were the gases, chiefly carbon monoxide and nitrogen, 
dissolved in the glass of the discharge vessel. As will be seen later, getting rid of 
these gases presented no difficulty. 

In Figure 1 the apparatus is shown with the barometric cut-offs and throttles 
set for the passage of the gas through the main circulation system consisting of 
the discharge vessel (10cm. diameter), the three-stage mercury vapour pump, 
. P,, and the reservoir (10cm. diameter). In L a high-tension arc could be 
maintained between a nickel anode and a molten lithium cathode for the puri- 
fication of the rare gas admitted to the system. By a suitable setting of the 
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appropriate cut-atts, the vessel K and the phosphoric oxide tube could be inserted 
into the circulation system between the reservoir and discharge vessels for the 
purpose of the quantitative analysis of steam-hydrogen-oxygen mixtures. Thus, 
with the platinum wire coil in K cold, the steam was removed with a corresponding 
change in the manometer level. Next, with the coil heated, hydrogen and oxygen 
in their combining proportions were removed, with a further change in the mano- 
meter level, the remaining gas being excess hydrogen, oxygen or known added 
inert diluent such as, for example, neon. The gaseous contents of the discharge 
vessel-reservoir circulation system could be well dried by circulation through the 
interposed phosphoric oxide tube, though that degree of rigid drying which 
exerts a marked effect in increasing the resistance of hydrogen to ignition could 
not be achieved. It is well known to workers in this field that, despite all pre- 
cautions, the use of mercury seals and cut-offs invariably introduces traces of 
moisture. 

The normal working pressure of gas in the reservoir during circulation was 
between about 7 and 10 mm., whilst that in the discharge vessel was maintained at 
between 0-026 and 0:03 mm. Hg by means of the pump, P,, and the appropriate 
throttle. The rate of circulation through the discharge vessel was rather over 
200 cm? per second ; this was measured in a manner which will be obvious from 
Figure 1. Owing to the fall in pressure at the throttle, the partial pressure of 
mercury vapour in the discharge vessel could not exceed 3 x 10-mm., the rate 
of circulation being more than sufficient to prevent diffusion of mercury from the 
McLeod gauge, cut-offs or pump towards the discharge vessel. 

Short sealed-off tubes were fused, as shown, into the discharge vessel and the 
1-2cm. bore tube, nearly 1 metre long, between the discharge vessel and the 
McLeod gauge entry tube. ‘These tubes served for the sealing in of probes 
projecting into the discharge or gas stream. ‘Three probe tubes, two of which 
are shown in Figure 1, were situated at 5, 20 and 50cm. respectively beyond the 
discharge vessel outlet. In all other respects Figure 1 is self-explanatory. 

The discharge was maintained by means of a high-frequency current in a 
solenoid, 13cm. in length and internal diameter, mounted coaxially with the 
discharge vessel, and consisting of seven turns of 4mm. diameter copper tubing. 
‘The current was obtained from a valve (Mullard S.W.9) oscillator in a modified 
Hartley circuit. ‘The valve filament was run at its rated value and, unless otherwise 
stated, the input to the oscillator was 160 watts. ‘The frequency was 26 Mc/s. 
as measured by a Southworth-Lecher (Southworth 1920) wire system. Further, 
the gas pressure in the discharge vessel was maintained at between about 0-026 
and 0-03mm. Hg, preferably nearer the former, the McLeod gauge readings 
being corrected for the experimentally determined equilibrium steam contents 
when moist gaseous mixtures were in use. 

Under these conditions the discharge in steam consisted of a bright rose-pink 
annulus coaxial with the discharge vessel, from the cool walls of which it was 
separated by a dark sheath about 2 or 3 mm. thick. In cross section the diameter 
of the bright core of the annulus was about 4 cm., the intensity of the glow decreas- 
ing considerably near the centre of the discharge tube. It will be evident that this 
| was aring discharge, the theory of which was outlined in the 13th Guthrie Lecture 
) by J. J. Thomson (1928). In such a discharge the currents form closed circuits 
| in the gas, though it is possible that, under the conditions of these experiments, 
the actual initiation of the discharge was effected by the electrostatic field across 
32-2 
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the ends of the solenoid (Mackinnon 1929, Esclangen 1934). It will be appreciated. 
that conditions in this ring discharge closely resemble those in the positive 
column. ; 

The evacuated discharge vessel was finally cleansed and gases absorbed in its 
walls sufficiently removed, first by prolonged baking at about 350° c. and then by 
maintaining the discharge in a stream of steam drawn from the vessel containing 
CO,-free water, past the appropriate throttle, through the cool discharge vessel 
and the pump P, backed by the ‘Hyvac’ pump. ‘This also served to free the 
water in its storage vessel from dissolved air. After evacuation following upon 
such cleansing, the main circulation route was restored and the system filled 
with the appropriate gas co the desired pressures in the discharge and reservoir 
vessels. 

A mixture of hydrogen and oxygen in combining proportions was frequently 
used in the experiments to be described below. It was. conveniently prepared 
in a high state of purity by drawing pure steam through the discharge and drying 
the products accumulated in the reservoir by subsequent circulation through the 
phosphoric oxide tube. Oxygen was obtained by heating recrystallized potassium 
permanganate and then washing with a caustic potash solution. Hydrogen was 
prepared by electrolysis of a saturated recrystallized baryta solution and, when. 
necessary, was subsequently freed from oxygen and hydrogen oxides by circulation 
past the heated platinum spiral K and through the phosphoric oxide tube. 


§3. EXPERIMENTAL RESULTS AND CONCLUSIONS 


Series 1—Dry and moist hydrogen. A puzzling feature of our work on the 
cathodic combustion of hydrogen was that, as has already been remarked, no 
evidence was found to justify the assigning of a direct réle in combustion to the 
atomic hydrogen, which was known to be produced in considerable quantities. 
Accordingly it was decided to examine first the behaviour of hydrogen, dried 
as well as possible by circulation over phosphoric oxide, in the electrodeless 
discharge. 

After a preliminary drying of the discharge tube walls by means of a prolonged 
circulation of dry hydrogen with the phosphoric oxide tube in the circuit, the 
filament K glowing and the discharge in action, the apparatus was evacuated and 
charged with dried hydrogen. It was then found that the hydrogen could be 
circulated through the discharge for long periods without any change-in pressure 
and that, in the early stages of circulation, the spectrum (Figure 2*) of the bluish- 
white discharge consisted mainly of the secondary spectrum due to neutral 
molecular hydrogen. ‘The hydroxyl bands (2808 and 30644.) and the more 
intense lines (H, 6563, Hg 4861 and H,, 4350.) of the Balmer series, due to neutral 
atomic hydrogen, were visible but, with the exception of H,, were feeble compared 
with those of the secondary series. Also in evidence were a few mercury lines,. 
in particular 2536, 3650 and 4358a. With increasing time of circulation and 
drying, during which precautions were taken to avoid shaking or otherwise 
disturbing the mercury in the McLeod gauge and the cut-offs associated with the 
main circulation system, the degree of dryness of the gas evidently improved until, 
- except for the brilliant secondary series and the few mercury lines mentioned above,, 
only the 3064. hydroxyl band-head region and the red H, line still remained. 


* For Figures 2 to 6 see Plates. 
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visible. ‘Thus it seems likely that, had it been possible to achieve a rigid drying 
of the hydrogen, the Balmer series and hydroxyl bands would have been eliminated. 

These results confirm, in all essential respects, the experimental findings of 
Wood (1920, 1921, 1922 a, b), Copeland (1930), Poole (1937) and other workers. 
‘Wood concluded that atomic hydrogen was copiously formed in the discharge 
passing through dry hydrogen and attributed the virtual non-appearance of the 
Balmer series to a catalytic recombination of the atomic hydrogen to the molecular 
form at the dry walls of the vessel. Other workers have since accepted this 
hypothesis. 

For the three following reasons, however, I have not been able to subscribe to 
this view without submitting it to the test of further experiment. In the first 
place, Wood and others worked with gas pressures of the order of 1 mm. ; it seems 
improbable, therefore, that any considerable proportion of hydrogen atoms 
formed in the discharge could arrive at the walls of the vessel without being 
excited and radiating, and thus giving rise to the Balmer series. Next, the catalytic 
inertness of cool, clean dry glass is well known, particularly to workers in the 
field of gaseous reactions in closed vessels. ‘Thus Strutt (1911) found that the 
greenish glow of active nitrogen in a clean glass vessel persisted for several minutes 
after excitation, although it was extinguished almost instantaneously on introducing 
into the vessel a superfically oxidized copper wire. Finally, although the quantum 
theory indicates the possibility of the dissociation of the hydrogen molecule by 
electron impact, the necessary excitation energy would be 18ev. (Mott and 
Massey 1933), and it is doubtful if such high energies were available in the ring 
discharge maintained under the conditions set forth above. ‘The matter is one 
which was put to the test of crucial experiment as follows. 

The pointed end of a glass needle was wetted to a length of about 3 mm. with 
diluted waterglass, dusted with yellow mercuric oxide powder, dried in vacuo 
-and sealed into the discharge vessel probe port. The tip of the needle was about 
level with the contour of the discharge vessel so that, if any hydrogen atoms were 
formed in the discharge, some would reach the mercuric oxide-coated needle tip 
directly. The circulation system with the phosphoric oxide tube in circuit and 
ithe wire K at red heat was then filled with hydrogen which was circulated for 
26 hours. ‘The discharge was then started in the, by this time, relatively well- 
dried, oxygen-free hydrogen, and the circulation continued. At first no change 
could be seen in the appearance of the mercuric oxide ; it was not until after 20 
minutes that the first signs of reduction, consisting in a blackening due to the 
formation of mercurous oxide, could be detected. ‘The drying tube was then by- 
passed out of the circulation system and a little oxygen admitted to the reservoir 
contents, whereupon, within a few seconds, the yellow changed to black oxide and 
this, in turn, was quickly further reduced to minute droplets of mercury which 
slowly evaporated. In another experiment a similar procedure was followed, 
except that instead of oxygen about 2 per cent by volume of steam was added to the 
reservoir vessel contents. Here again the same remarkable acceleration in the rate 
of reduction of the mercuric oxide occurred. As will be shown later, the agent 
responsible for the reduction of both mercury oxides was atomic hydrogen. 

It is difficult to resist drawing from these experimental results the conclusion 
that, apart from traces due to imperfect drying and revealed spectrographically, 
atomic hydrogen is not formed when an electrodeless discharge passes through 
dry hydrogen. Hence, the hitherto generally accepted view that the hydrogen 
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molecule is atomized by the discharge in dry hydrogen and that the absence of 
the Balmer spectrum is to be accounted for by a recombination to the molecular 
state at the surface of the discharge vessel is untenable. These experiments, 
of course, afford no information as to the relative activities of dry and moist glass 
surfaces in causing the recombination of hydrogen atoms to molecules ; but. 
they do show that dry hydrogen is not dissociated by the discharge, and that the 
minute amount of atomic hydrogen formed is attributable to that trace of moisture 
which can hardly be eliminated in an experiment of this kind carried out in an 
apparatus equipped with barometric mercury cut-offs. 

Series 2—The steam equilibrium. Preliminary experiments in which steam 
had been circulated through the discharge, primarily with a view to outgassing 
the discharge vessel, had shown that much of the steam was decomposed and that. 
the only other products received into the reservoir vessel consisted of hydrogen 
and oxygen in their combining proportions. In further such experiments steam 
was drawn through the discharge, and the products accumulated to a pressure of 
about 10mm. in the previously evacuated reservoir, the cut-off between this and 
the discharge vessel being closed. On interrupting the steam supply and restoring 
the main circuit, no further change in pressure occurred, even after several hours’ 
circulation. Analyses immediately before resumption of, or after, prolonged 
circulation revealed no change in the steam concentration, the remaining products 
always consisting of hydrogen and oxygen in their combining proportions. It is 
clear, therefore, that the steam dissociation equilibrium was attained in a single 
passage of the steam through the discharge and, further, that no gas was lost by, 
for example, combination with mercury. 

In the course of these and other exploratory experiments it was also found 
that the extent of this dissociation was much affected by the state of the surface 
of the walls of the discharge vessel. For example, with clean walls and starting 
with either steam or the equivalent mixture of hydrogen and oxygen, the products 
at equilibrium contained between 16-5 and 19 per cent by volume of steam. 
On repeating such an experiment, but with a fine (37s.w.G.) copper wire protruding 
into the discharge through the appropriate port, the wire glowed and the equili- 
brium steam concentration rose progressively in the course of five two-hour runs 
to between 38 and 40 per cent. By this time the walls of the discharge vessel. 
near the wire were faintly discoloured, presumably by evaporated copper which 
may or may not have been oxidized. The glowing of the wire was due to some 
energetic surface reaction, and not merely to eddy-current heating, as was shown 
by the fact that on reducing the pressure in the discharge vessel, by closing the 
throttle and thus extinguishing the gas discharge, the wire at once ceased to glow.. 
The wire was then removed and the port sealed off, whereupon further runs, in 
which equilibrium was approached sometimes from one side and sometimes from 
the other, continued to yield equilibrium steam concentrations approaching: 
40 per cent as before. Thus there was no tendency for this concentration to 
increase, as had been the case when the copper wire protruded into the discharge. 

The discharge vessel was now cut out of the apparatus, cleaned with fuming 
nitric acid followed by chromic acid, washed with distilled water and finally dried 
with warm filtered air. The vessel was then re-sealed into the circulation 
apparatus, using dry filtered air for glass-blowing purposes. It was then found. 

‘that this cleaning had restored the discharge vessel to its original state, equilibrium 
steam concentrations of about 18 per cent being once more obtained. 


Steam in the Ring Discharge 47L 


These results point to a promotion, by the copper stain on the discharge: 
vessel walls, of some surface reaction or reactions which, either directly or through 
the products formed, favoured the formation of steam. 

The reactants concerned must be relatively long-lived, because they survive 
passage from the discharge region to the vessel walls across the intervening dis- 
charge-free dark zone. Further, bearing in mind the evidence afforded by the 
glowing copper wire of the highly energetic nature of the reaction, and the testimony 
of the mercuric oxide probe experiments with moist hydrogen to the arrival of 
atomic hydrogen at the discharge vessel walls, it seems likely that the surface 
reaction is the recombination of atomic to molecular hydrogen, and that this may 
be, indeed, the dominating, reaction increasing the equilibrium steam concen- 
tration. In other words, the formation of atomic hydrogen seems, somehow or 
other, to hinder the formation of steam. 

Further light was thrown on these matters by the spectrograms (Figure 3): 
ot the discharge in steam and its decomposition products. These are striking in 
that, apart from the relatively weak mercury spectrum, they consist of the Balmer 
series (3771, 3797, 3835, 3890, 3970, 4102, 4340, 4861 and 6562a.), the hydroxyl 
bands with heads at 2808 and 3064 4., and the rather weak lines, 3947 and 4368 a., 
due to atomic oxygen. ‘The secondary molecular hydrogen spectrum is absent. 
Further, the presence or otherwise of copper deposited on the discharge vessel 
walls had no appreciable effect on the appearance of the spectrograms or on the 
relative intensities of the hydroxyl band and Balmer spectra, despite the profound. 
effect exerted by the copper stain in the equilibrium steam concentration. 

Wood and other workers have obtained spectra under conditions comparable 
to those outlined above, except that, instead of pure steam, they used moist: 
hydrogen, the steam concentration being in general about 24 per cent. Apart 
from the absence of all atomic oxygen lines, their spectra are essentially similar 
to those obtained in the present experiments. Thus, while Wood and others. 
have shown that copious supplies of long-lived atomic hydrogen are obtained. 
from the discharge in moist hydrogen, the above experiments demonstrate that 
atomic hydrogen is also freely produced when steam alone is fed to the discharge. 

Series 3—Steam-hydrogen mixtures. In order to gain further insight into. 
the mechanism of the formation of atomic hydrogen, experiments were now 
carried out in which a wide range of mixtures of hydrogen with steam, or with 
hydrogen and oxygen in their combining proportions, were circulated through 
the discharge. Before each such experiment, however, the circulation system. 
was prepared in the following manner. Glass needles, the points of which had 
been coated with ash (mainly tin oxide and lime) by calcining a thin smear of 
sealing wax, were sealed into the three ports with their tips projecting into the 
tube between the discharge vessel and the McLeod gauge. Another ash-coated 
probe was sealed into the discharge tube port with its tip level with the contour 
of the vessel. Dry oxygen was then circulated under standard pressure conditions. 
through the discharge. The reaction product reacted with mercury vapour to 
deposit in a few minutes broad rings or, rather, bands of yellow mercuric oxide 
about the openings leading to the McLeod gauge and the rare gas purifier, and 
just above the entry to the pump, P,. The apparatus was then evacuated and the 
circulation system, with the phosphoric oxide tube in circuit, swept out and finally 
charged with roughly dried hydrogen, which was then still further dried by pro- 
longed circulation during which the platinum spiral, K, was glowing. ‘Throughout 
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these operations the mercuric oxide films remained unchanged. ~The discharge 
was now started and spectrographed. Apart from the strong secondary molecular 
spectrum and the usual lines due to traces of mercury, only the H, (moderate to 
weak), H, (weak) and H, (very weak) lines of the Balmer series and the hydroxyl 
3064 a. band-head region were visible. The circulating gas had no visible effect 
on any of the probes coated with sealing-wax ash ; likewise, the mercuric oxide 
films near the McLeod gauge and the pump were unchanged, even after two hours’ 
circulation. 

The apparatus was now evacuated, and the normal circulation route restored 
and filled to the usual pressures with the required mixture of hydrogen and steam. 
After a sufficiently prolonged circulation to ensure thorough mixing of the gases, 
the discharge was started. The range of mixtures examined extended from 
steam alone to steam (or the equivalent hydrogen-oxygen mixture) diluted with 
up to 95 per cent of hydrogen. Spectrograms were obtained in each case under 
uniform conditions of exposure and development. ‘The sealing-wax ash probe 
inserted through the discharge vessel port was withdrawn after the first experiment 
and was not replaced. 

The results of this series are illuminating. In the first experiment, carried 
out with steam alone, the mercuric oxide films, including that near the mercury 
vapour pump inlet, about 14 metres beyond the discharge vessel, were completely 
reduced within a few minutes of starting the discharge. ‘The progress of such 
reduction could, however, be immediately arrested by interrupting the discharge. 
Further, the particles of ash on the discharge vessel probe and on the two probes 
nearest to the discharge vessel flashed into incandescence immediately the 
discharge was started, but ceased to glow when this was interrupted. From time 
to time glowing fragments of ash would fly off the probes, particularly that in the 
discharge vessel, and these continued to incandesce after lodging on the wall of 
the vessel or in the tube beyond it. Indeed, owing to this stripping-off of the ash, 
the section of tubing carrying the first two ports beyond the discharge vessel had 
to be removed for cleaning after each experiment and fresh ash-coated probes 
inserted in preparation for the next. As stated above, the discharge vessel 
probe was removed and not replaced after the first experiment. 

The incandescence of the ash on the second probe beyond the discharge 
vessel was less brilliant than that of the first probe and, in the case of the third 
probe some 50cm. beyond the discharge vessel, an ash particle would sometimes 
glow, but only in mixtures containing initially between about 30 and 60 per cent 
of hydrogen. 

The spectrograms were similar to those obtained with pure steam (Figure 3), 
except that the atomic oxygen lines, 3947 and 4368a., seemed to weaken with 
increasing hydrogen concentration ; they no longer appeared in a spectrogram 
from the discharge in a mixture of steam and hydrogen in equal volumes. Visual 
comparison of the spectrograms, more particularly of the hydroxyl band 2808 a. 
with H3797 and H3835a., and of the hydroxyl band 3064 a. with H3890 and 
H3970a., revealed no change with increasing hydrogen concentration in the 
relative intensities of the hydroxyl band and Balmer series spectra. 

Finally, it. was found that the steam equilibrium constant decreased with 
increasing dilution of the mixture by hydrogen. For example, with steam alone 
the constant K =(H,O)?/[(H,)?. (O2)] =0-46, but was 0-38 and 0-25 in the case of 
mixtures diluted with 55 and 69 per cent of hydrogen respectively. Thus, 
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although dilution with hydrogen did, in fact, greatly reduce the dissociation of 
‘the steam by the discharge, it did not do so to the extent required by the law of 
mass action. 

The reason is not far to seek. The results of these experiments as a whole 
now afford overwhelming evidence that the first step towards equilibrium in 

.steam or in a steam—hydrogen mixture is the direct dissociation of the steam into 
hydroxyl and atomic hydrogen. The molecular hydrogen as such is not dis- 
-sociated by the discharge and, when steam is present in adequate supply (which 
may possibly be considerably less than 5 per cent), either the hydrogen is not 
-excited or it reacts in some way before it can radiate. This dissociation of steam 
must be practically irreversible in the gaseous phase, because quantum energy 
‘considerations forbid the reverse reaction by which atomic hydrogen and hydroxyl 
might recombine to form steam, except as a result of necessarily infrequent three- 
body collisions. Similar considerations apply in the case of the union of hydrogen 
atoms to the molecular state. On the other hand, the life of the hydroxyl radical 
‘(Bonhoeffer and Pearson 1931, Geib and Harteck 1934, Oldenberg 1937) can be 
terminated by two-body collision processes such as, for example, 20H =H,O +O, 
and OH+H,=H,0+H. ‘Thus the fall in the steam dissociation equilibrium 
constant with increasing dilution by hydrogen is to be accounted for by the 
practically irreversible formation, and hence escape, of atomic hydrogen from the 
reaction zone. 

There is evidence from the above experiments that the hydroxyl reacted in 
the discharge in both the ways suggested. ‘That the first reaction occurred is 
supported by the appearance of oxygen atom lines in the spectrograms from the 
discharge in steam, and by their weakening and final disappearance on progressive 
dilution with hydrogen. As for the second reaction, in which hydroxy] interacts 
with hydrogen, the likelihood of its occurrence is supported by the following 
facts and considerations. In the first place, the relative intensities of the hydroxyl 
and Balmer spectra (Figure 3) were independent of hydrogen dilution of steam 
over a wide range. ‘This strongly supports the previous conclusion, that the sole 

‘primary source of atomic hydrogen is the direct dissociation of steam in the 
discharge. Next, the fact that the third sealing-wax ash probe glowed only in 
mixtures of steam diluted by a roughly equal proportion of hydrogen shows that 
the yield of atomic hydrogen was increased by such dilution. Hence the diluent 
hydrogen must itself contribute in some way to the ultimate yield of atomic 
hydrogen. Since, as has already been shown, atomic hydrogen cannot have been 
produced by direct dissociation of molecular hydrogen, the probable source of 
this contribution lies in the interaction in the discharge of hydroxyl and molecular 
hydrogen, whereby atomic hydrogen and steam are formed. ‘T’his view received 

strong support from the results, which showed that the recombination of atomic 
hydrogen at the copper-stained walls of the discharge vessel increased the steam 
equilibrium concentration by returning molecular hydrogen to the discharge and 

‘thus reducing the amount escaping in the atomic state from the reaction zone. 
The possibility of a surface reaction between hydroxyl and atomic hydrogen in 

_ similarly affecting the steam equilibrium will be considered later. 

The results of this series of experiments furnish ample evidence of the signi- 
ficance to be attached to the probe and mercuric oxide reactions. Thus, when 

‘the discharge was fed with steam, steam diluted with hydrogen, or the equivalent 
hydrogen—oxygen mixtures, the mercuric oxide films were rapidly reduced, and 
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the sealing-wax ash probes incandesced with a brilliancy which increased as the- 
discharge vessel was approached. Further, all these phenomena ceased imme- 
diately the discharge was interrupted, irrespective of whether the circulating gas 
was steam (or the equivalent hydrogen-oxygen mixture) alone or diluted with 
hydrogen. Therefore the agent reducing the mercuric oxide could not have been 
molecular hydrogen, nor could the reaction catalysed by the ash have been the 
union of hydrogen and oxygen to steam. In addition to the normal equilibrium 
products entering the reservoir vessel, the only other products for which there is 
evidence of formation in the discharge are atomic hydrogen, hydroxyl and hydrogen 
peroxide (Finch and Mahler 1931). Atomic oxygen can be excluded because, 
although revealed spectrographically in the case of steam traversing the discharge, 
it was eliminated by sufficient dilution of the reaction mixture with hydrogen.. 
Hence, whereas the reduction of the mercuric oxide films can only be attributed to 
atomic hydrogen, the incandescence of the ash-coated probes might be ascribed 
either to the formation of molecular hydrogen or to the interaction of atomic 
hydrogen with hydroxyl to form steam. The heat evolutions of these reactions 
are similar, being about 103,000 kcal/gm. mol (see Steacie 1946). Whilst there 
is no doubt about the longevity of atomic hydrogen, the situation in regard to the 
life of hydroxyl, though known to be much shorter than that of atomic hydrogen,. 
is less clear. Thus, while Bonhoeffer and Pearson (1931) were unable to find the 
radical in the gas flowing from a discharge through water vapour, Frost and 
Oldenberg (1936), using a 21-foot grating, detected hydroxyl for as long as 
4 second after their discharge in steam was interrupted. Thermal conditions may, 
however, have contributed to this result. The reduction of the mercuric oxide 
probe in the discharge vessel port, when the discharge took place in moist hydrogen, 
would be difficult to explain if hydroxyl arrived at the probe, and hence also at 
the surface of the vessel, at a rate comparable with that of the atomic hydrogen.. 
Thus the evidence so far supports the view that the incandescence of the ash-coated 
probes was due to the recombination of atomic hydrogen. This issue was further 
tested in the next series of experiments. 

Of the three ash-coated probes, all of which were much nearer to the discharge 
vessel than were the mercuric oxide films, the first two incandesced regularly, 
whilst the third did so spasmodically, and then only under the most favourable- 
conditions. It would seem that, whilst the reduction of the mercuric oxide film 
formed a sensitive test for the presence of atomic hydrogen, the sealing-wax ash 
probes could respond only to a relatively high concentration. 

Series 4—Further probe experiments. In the first experiment of this series. 
a closely wound spiral nearly 1cm. in diameter, of 28 s.w.G. copper wire, was 
inserted into the exhaust tube leading from the discharge vessel. The spiral, 
previously cleaned and reduced by hot hydrogen, extended from near the second. 
port, which was fitted with a mercuric oxide-coated probe, to within 0-5 cm. 
trom the discharge vessel. ‘The equilibrium steam concentration was then 
determined, starting with steam or the equivalent hydrogen—oxygen mixture. 
The mean steam yield from eight such runs was 19-2 per cent, a value which did 
not differ significantly from that (18-9 per cent) afforded by a similar series of runs. 
carried out before inserting the spiral. The presence of the copper spiral, 
however, practically inhibited the reduction of the mercuric oxide, thus showing 
that the bulk of the atomic hydrogen leaving the discharge vessel was, on passing: 
through the spiral, transformed to the molecular state. This, together with the 
fact that the steam yield was virtually unaffected by the presence or otherwise of 
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the spiral, proves conclusively that hydroxyl could not have survived long enough 
after leaving the discharge to reach the spiral in appreciable quantities ; otherwise 
the equilibrium steam concentration must have increased, which it did not do. 
Hence these results confirm that the ash incandescence observed in the previous. 
experiments was due solely to recombination of atomic hydrogen. 

The spiral was now removed and the corresponding tube section cleaned and 
re-fused into the apparatus. ‘The glass pinch of a tungsten filament lamp was. 
then sealed into the first probe port beyond the discharge vessel, the lamp filament 
having been replaced by a 2:5cm. long fine (37s.w.c.) platinum wire and the 
current-carrying posts shrouded with quilltubing. The wire projected in the form 
of a hoop into the path of the gas stream from the discharge vessel. 

In the next experiment a mixture of hydrogen and oxygen was circulated under 
normal conditions but without the discharge, the phosphoric oxide tube being 
included in the circulation path. After drying the gas in this manner until the 
manometer readings were steady, the platinum filament was brought to a bright 
red heat and circulation continued. The McLeod gauge readings then fell from 
the steady value of 0-027 mm., obtained before heating the filament, to 0-019 mm., 
thus clearly indicating that steam was being formed as a result of a surface union 
of hydrogen and oxygen on the heated platinum wire. This was further confirmed 
by a progressive fall in the reservoir manometer readings as circulation proceeded. 
On interrupting the heating current, the wire immediately cooled and the McLeod 
gauge reading rose to its previous steady value of 0:027mm._ Evidently the rate 
of heat evolution during the formation of steam at this low pressure was not 
sufficient to compensate for the heat losses from the wire, which in consequence 
cooled and ceased to catalyse the reaction. Immediately the discharge was 
started, however, the platinum filament began to glow brightly, but cooled at 
once when the discharge was switched off. ‘That the heating of the wire must 
have been due solely to the recombination of atomic hydrogen, and not to heating 
of the filament by the high-frequency field exciting the discharge, was shown by 
tha fact that the wire cooled immediately the discharge was interrupted by closing 
the throttle admitting gas into the discharge vessel. 

The difference in behaviour of the platinum wire in atomic hydrogen and in 
a mixture of hydrogen and oxygen in equivalent proportions finds a natural 
explanation in the heat evolution of the appropriate reactions, the heat of forma- 
tion of molecular hydrogen being 103,000 kcal/gm. mol, while that of steam is 
only about 59,000. 

Further experiments were carried out with, in place of the platinum filament, 
a probe consisting of a straight length of either 36 s.w.c. tungsten or 37 s.w.c. 
copper wire protruding into the gas stream from a glass rod sealed into the port 
adjacent to the discharge vessel. When steam was fed to the discharge, these 
wires glowed dull red. With a mixture of steam and hydrogen in equal pro- 
portions, the copper wire became bright red and the tungsten wire was even hotter. 

Wood (1922a) also found that a tungsten wire glowed in atomic hydrogen ; 
he did not however obtain such an effect with a platinum wire. ‘This may have 
been due to the state of the surface of his wire. It is a common experience that 
catalysing the combustion of a mixture of hydrogen and oxygen at a bright red 
heat not only cleans and activates a platinum surface but also roughens it. 

Series 5—Steam-—mercury vapour mixtures. At the conclusion of the 
experiments of Series 3 the question was left open as to whether or not the effect 
of the copper-staining of the discharge vessel walls in promoting a high yield of 
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steam might be due, not only to the returning of molecular hydrogen to the reaction 
zone, but also to a surface reaction between atomic hydrogen and hydroxyl. In 
previous work on cathodic combustion (Finch 1935) it had been suspected that 
mercury vapour facilitated the combustion of hydrogen. ‘This, and the fact that 
mercury vapour can play an important réle in photosynthesized reactions involving 
hydrogen (Taylor and Bates 1926, Rieke 1936, Olsen 1938), suggested that its 
use might throw light on the present problem. Accordingly arrangements were — 
made to examine the effect of mercury vapour additions to steam passing through 
the ring discharge. 

About 0:5cm* of mercury was distilled into the discharge vessel port tube 
through a branch tube, which was then sealed off close to the port tube. A few 
turns of resistance wire wrapped round the port tube and immersed in an oil- 
filled jacket enabled the mercury to be warmed. An ash-coated probe was inserted 
through the first port beyond the discharge vessel. Mercuric oxide films were then 
formed, in the manner previously described, at the usual places where the circu- 
lating gas stream came in contact with stagnant mercury vapour. A similar 
‘deposit also formed about the entry to the discharge-vessel port tube. 

Steam (or the equivalent hydrogen—oxygen mixture) was then circulated along 
the usual path under the normal pressure distribution of 0-026 mm. in the discharge 
vessel and 7 or 10 mm. in the reservoir, according to the nature of the gas mixture 
used. ‘The mercury in the discharge port vessel was then gently heated and 
maintained at between 40 and 45°c. When steady conditions had been attained, 
the discharge was started. Instead of a rose-pink ring, it now formed an apparently 
uniform bluish-green sphere. A rough estimate afforded by the increased damping 
of the oscillation circuit suggested that the partial pressure of mercury vapour 
was between 0-005 and 0-012mm., probably nearer the former. Despite this, 
no condensation occurred anywhere between the discharge vessel and the pump 
though, with time, a barely perceptible brownish band formed on the discharge 
vessel walls coaxial with the solenoid. 

The mercury vapour exerted a profound effect on the results. The ash- 
coated probe did not incandesce and, although the mercuric oxide film near the 
discharge vessel port had disappeared within five minutes of starting the discharge, 
the more remote films were only very slowly reduced. Also, the equilibrium 
steam concentration, which previously had been between 18 and 19-4 per cent, 
was now practically 100 per cent, irrespective of the direction of approach to 
equilibrium. Finally, the spectrogram (Figure 4) of the discharge differs 
strikingly from that obtained with steam or steam—hydrogen mixtures (Figure 3). 
Apart from the strong, well-developed mercury line spectrum, its main features 
are the 2808 and 30644. hydroxyl bands, the latter being intense, and the three 
mercurous hydride line groups 3500, 4017 and 4219. Neither the secondary 
hydrogen nor the Balmer spectra appear. | 

These results show that the effect of the mercury vapour was largely to 
suppress the formation of atomic hydrogen, thereby shifting the equilibrium so 
far towards the steam side that the reaction products consisted almost wholly of 
steam. ‘hus the evidence suggests that steam is reduced, reversibly, by mercury 
to mercurous hydride and hydroxyl. In view of the high intensity of the 2536. 
mercury line [Hg(6°P,)] the reduction of mercury by hydrogen to mercurous 
hydride and atomic hydrogen is also likely to occur. Again, since mercury is 
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oxidized by ozone to mercuric oxide, and this in turn is reduced in the cold by 
atomic hydrogen and by molecular hydrogen at about 120°c. (Manchot and 
Kampschulte 1907, 1910), further probable reactions are the oxidation of mercurous 
hydride to mercuric oxide and hydroxyl, and the reduction of mercuric oxide by 
atomic or molecular hydrogen to mercury and hydroxyl or steam respectively. 
In addition, the short life of hydroxyl is terminated in the manner set forth in 
Series 3 and 4 of the above experiments. 

The question is still left open as to whether the survival life of hydroxyl is 
sufficient for it to interact with atomic hydrogen at the walls of the discharge 
vessel ; the results of Series 4 show, however, that this reaction does not occur 
at surfaces beyond the discharge vessel. 

Series 6—Dry and moist oxygen mixtures. The behaviour of oxygen, both 
dry and moist (or diluted with hydrogen), in various types of discharges, including 
the ring discharge, has been the subject of much experimental work. In general, 
a broad measure of agreement has been reached as to findings and conclusions. 
In one respect, however, a conflict of views does exist. Harteck and Kopsch 
(1931), Geib and Harteck (1934), Geib (1938) and others incline to the belief 
that, like atomic hydrogen, atomic oxygen has a relatively long life. On the other 
hand, Strutt (1910) and Finch and Bradford (1934) froze out ozone from the gas 
issuing from a discharge fed with dry oxygen, but found no evidence of the survival 
of atomic oxygen. All these workers appear to have used discharges between 
electrodes. Copeland (1930), however, used an electrodeless discharge, and 
concluded that atomic oxygen not only has a life comparable with that of atomic 
hydrogen but is also formed in the discharge by a mechanism similar to that 
postulated by Wood, whereby a poisoning of the discharge vessel walls by moisture 
was supposed to be necessary for the survival of the atomic species. It has, however, . 
already been mentioned above that mercuric oxide films were formed after a 
brief circulation of oxygen through the discharge. ‘These films were as readily 
obtained with dry as with moist oxygen. Thus if, as Copeland held, they were 
due to the interaction of mercury vapour and atomic oxygen, their formation with 
dry oxygen conflicted with his hypothesis. It was therefore decided to examine 
the matter further. 

Owing to the difficulty of removing stains, caused by the experiments of 
Series 5, from the discharge vessel without etching the glass and thus impairing 
its smooth surface, a new vessel was fitted. It had previously been cleaned in 
the usual way but not baked out or degassed by circulating steam through the 
discharge. After fitting the first port beyond the discharge vessel with an ash- 
coated probe and the discharge vessel port with a 36s.w.c. tungsten wire protrud- 
ing into the discharge region, the normal circulation system was filled with oxygen 
previously washed with concentrated caustic potash and roughly dried by a single 
passage over phosphoric oxide. The gas was then circulated through the dis- 
charge under the normal conditions of pressure. 

Rapid formation of the mercuric oxide films again occurred and was followed 
more slowly by a “‘tailing”’ of the mercury in the pump boiler. Neither the ash- 
coated nor the tungsten probe showed any signs of incandescence. ‘The 
‘tailing’ and mercuric oxide films were easily removed by reduction, in the manner 
previously described, with atomic hydrogen formed from a hydrogen-steam 
mixture in equal proportions. A spectrogram, taken about half an hour after 


478 G.I. Finch 


starting the discharge, is reproduced in Figure 5. It reveals, inter alia, a fairly 
strong atomic oxygen spectrum, faint hydroxyl bands, the first three lines of the 
Balmer series, and prominent bands due to carbon monoxide. 

On repeating this experiment, a spectrogram similar to Figure 5 was obtained, 
but without the carbon monoxide spectrum. It would appear, therefore, that in 
the first experiment of this series the circulation of oxygen through the discharge 
had effectively out-gassed the discharge vessel walls. A fresh charge of oxygen 
was circulated for several hours past the heated spiral K and through the drying 
tube before starting the discharge, and a further spectrogram taken. ‘This was 
again similar to Figure 5, except that, in addition to the absence of the carbon 

_monoxide spectrum, the 3064, but not the 2808 a., hydroxyl band-head region 
was visible. H, of the Balmer spectrum was very faint. 

It will be evident at once that, in affording the primary (atomic) spectrum, 
the behaviour of dry oxygen in the ring discharge differed radically from that of 
dry hydrogen which, in the main, yielded the secondary spectrum, with weak 
hydroxyl bands and a few faint lines of the Balmer series attributable to traces 
of moisture. 

It is noteworthy that, despite the spectrographic evidence of the production 
of atomic oxygen in the discharge, the tungsten and ash-coated probes remained 
cool throughout. Since the heat of formation of molecular oxygen (>117,000 
kcal/gm. mol.) is greater than that of hydrogen, it is plain that little or no atomic 
‘oxygen survived to reach the ash probe, only 5cm. beyond the discharge vessel, 
or even the tungsten wire probe. It is most unlikely, therefore, that the mercuric 
oxide films, which had formed so rapidly, could have been due to the interaction 
of mercury vapour and atomic oxygen, as Copeland has supposed. Rather must 
this phenomenon be ascribed to some other active modification of oxygen. 

In the next experiments the effect of the dilution of oxygen by steam was 
-examined. A wide range of mixtures was used containing between 5 and 96 per 
cent steam, the remainder being oxygen. ‘The results have so much in common 
that it will suffice to summarize them as follows: 

(i) Mercuric oxide films were formed, as with dry oxygen, at the usual sites, 
-and ‘‘tailing”’ in the pump also occurred. The rate of formation of the mercuric 
oxide, as estimated from the growth in density of the film near the McLeod gauge, 
was roughly proportional to the oxygen concentration. 

(ii) The ash-coated probe incandesced to some extent in all cases, but did so 
most vigorously in mixtures rich in steam. ‘The tungsten wire probe glowed 
‘brightly in all cases. 

(iii) The spectrograms reveal a progressive weakening of the atomic oxygen 
‘spectrum with increasing steam concentration, but a few lines, notably 3947 
-and 43684., remain, even in the case of pure steam. ‘The hydroxyl band and 
Balmer spectra are prominent throughout. 

It has already been shown in the case of dry oxygen that the formation of the 
mercuric oxide films was due to some active form of oxygen, other than atomic 
oxygen, which was incapable of exciting either the ash-coated or tungsten wire 
‘probe to incandescence. Hence, the fact that this oxidation process was slowed 
down by progressive dilution of the oxygen with steam strongly suggests that in 
‘the cases of both dry and moist oxygen the same agent was responsible for the 
oxidation of mercury. The results obtained by Strutt (1911) and Finch and 
Bradford (1934), taken into consideration with the fact that ozone is known to form 
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yellow mercuric oxide with mercury (Manchot and Kampschulte 1907, 1910), 
justify the assumption that ozone was the long-lived, active form of oxygen 
produced in these experiments. 

Since the ash probe and, toa lesser extent, the tungsten wire incandesced more 
brightly with increasing dilution with steam, it is evident that this could only have 
been due to recombination of atomic hydrogen. That the phenomena of the 
probe incandescence and mercuric oxide formation occurred simultaneously, 
even in mixtures containing as little as 5 per cent of oxygen, suggests that atomic 
hydrogen and ozone interact in the gas phase only relatively slowly, if at all, under 
the conditions of our experiments. 

In the light of these results we are forced to conclude that, unlike atomic 
hydrogen, atomic oxygen does not survive as such on escaping from the discharge. 
‘There is a ready explanation for this. ‘Thomson and Thomson (1928) point out 
that, unlike hydrogen, the oxygen molecule readily captures an electron to form 
a negative ion. It is to be expected that the oxygen atom is even more efficient 
in this respect. ‘There is ample opportunity for such capture to occur within the 
discharge, where the electron population is of the same order as that of the mole- 

cules present (Thomson 1928). Under these circumstances, two-body collision 
reactions, such as the recombination of oxygen atoms to molecular oxygen, or 
of oxygen atoms with molecules to form ozone, are no longer subject to quantum- 
mechanical restrictions, provided one or both of the reactants is a negative ion. 

Series 7—Neon-steam mixtures. ‘This series of experiments forms, in a sense, 
a continuation of Series 6, the oxygen component being, however, replaced by a 
rare gas. It was hoped that a comparison of the results of the two sets of experi- 
ments would throw further light on, zuter alia, the apparent strange reluctance of 
atomic hydrogen and ozone to interact. The experiments were carried out with 
mixtures covering a range from well-dried neon to pure steam. 

The essential results can be Summarized briefly by stating that (i) the intensity 
of the incandescence of the ash and tungsten wire probes increased with increasing 
steam concentration in a manner which closely paralleled that observed in the 
case of the corresponding oxygen-steam mixtures and, (ii) unlike the case of 
steam-diluted oxygen, the spectrograms (Figure 6) show that on increasing the 
steam concentration the neon spectrum was quickly suppressed. For instance, 
in a mixture of neon and steam in equal proportions only the hydroxyl and Balmer 
series spectra are visible. 

The parallel behaviour of the probes in the two series supports the conclusion 
that in the case of both oxygen and neon-steam mixtures incandescence of the 
probes was caused solely by atomic hydrogen and, further, that atomic hydrogen 
and ozone are, for some unknown reason, slow to react with each other. It is 
also not difficult to see that, except for the formation or otherwise of mercuric 
oxide, the experimental results which led Copeland to suppose that atomic oxygen 
of long life is formed by the action of the electrodeless discharge on moist oxygen 
could equally well have been obtained with moist neon. 


§4. SUMMARIZED FINDINGS AND CONCLUSIONS 


The behaviour of steam or a mixture of hydrogen and oxygen in combining 
| proportions, and their mixtures with hydrogen or oxygen, ranging from dry 
__ hydrogen to dry oxygen, have been studied in the ring discharge. The effect of 
| mercury vapour and of neon on the steam equilibrium was also examined. 
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The principal facts established were as follows : 

(i) Dry hydrogen. ‘The secondary hydrogen spectrum was strongly emitted.. 
The hydroxyl-band and Balmer spectra were weak and tended to disappear with 
progressive drying. Apart from traces attributable to incomplete drying, no 
agent capable of reducing yellow mercuric oxide at room temperature reached 
the walls of the discharge vessel. 

(ii) Moist hydrogen. On adding to dry hydrogen about 1 per cent of oxygen 
or 2 per cent of steam (and possibly much less) the spectrum changed radically, 
the secondary hydrogen spectrum being replaced by the hydroxyl-band and 
Balmer spectra. Mercuric oxide at the wall of the discharge vessel was rapidly 
reduced. 

(iii) Steam. 'The emission consisted mainly of the hydroxyl-band and Balmer 
spectra, together with a few lines due to atomic oxygen. 

The stable end-products obtained from steam (or the equivalent hydrogen— 
oxygen mixture) consisted of steam, hydrogen and oxygen, the last two in their 
combining proportions. With a clean discharge vessel, the final steam concen- 
tration was between 16-5 and 19 per cent, but partial contamination of the 
vessel’s walls by a barely visible film of copper sufficed to more than double the 
steam yield, which, however, returned to its previous value after removal of the 
copper film. 

Glass needles, coated with sealing-wax ash (SnO and CaO) and inserted into 
the exhaust tube through ports situated 5 and 20 cm. beyond the discharge vessel, 
incandesced brightly. Films of yellow mercuric oxide deposited in the exhaust 
tube were rapidly reduced to metallic mercury. Copper, tungsten or platinum. 
wire probes in the discharge vessel or in the exhaust tube glowed. The insertion. 
of a copper wire spiral into the exhaust tube between the discharge vessel and an: 
ash-coated probe suppressed the incandescence of the probe and practically 
stopped the reduction of the mercuric oxide films. The spiral did not appreciably 
affect the equilibrium steam concentration. With the ring discharge interrupted, 
an electrically heated platinum wire probe, situated in the exhaust tube near the 
discharge vessel, catalysed the formation of steam from the equivalent hydrogen— 
oxygen mixture, but ceased to do so and cooled when the heating current was 
stopped. The wire glowed spontaneously, however, immediately the ring dis-- 
charge was restored. 

(iv) Steam-hydrogen and the equivalent hydrogen—oxygen mixtures. 'The 
hydroxyl-band and Balmer spectra were prominent in all mixtures contaiiing up: 
to 95 per cent (and possibly more) diluent hydrogen, and their relative intensities 
appeared to be independent of such dilution. On increasing dilution with hydro- 
gen, the atomic oxygen lines observed with steam alone weakened and eventually 
disappeared. 

Dilution with hydrogen reduced the dissociation of steam, but not to the extent _ 
required by the law of mass action. 

The ash-coated probes incandesced and the mercuric oxide films were reduced 
in all mixtures examined, incandescence being most pronounced in the case of 
mixtures containing hydrogen and steam in roughly equal proportions. 

(v) Steam—mercury vapour. Apart from an intense mercury spectrum, the 
emission from the ring discharge in a mixture of steam (or the equivalent hydrogen— 
oxygen mixture) with roughly one-third its volume of mercury vapour consisted 
chiefly of the hydroxyl-band and mercurous-hydride spectra. The Balmer and 
secondary hydrogen spectra were absent. The ash probes did not incandesce.. 
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The mercuric oxide films in the exhaust tube were only slowly reduced. In 
contrast to the steam yield of about 18 per cent obtained with pure steam, the 
equilibrium steam concentration rose to nearly 100 per cent. 

(vi) Dry oxygen. The emission consisted mainly of the atomic oxygen 
spectrum, with faint hydroxyl-bands and a few faint lines of the Balmer series 
attributable to traces of moisture. 

Yellow mercuric oxide was formed in the exhaust tube at all those sites where 
the gas stream issuing from the discharge came into contact with stagnant mercury 
vapour. ‘The wire and ash-coated probes remained cool. 

(vii) Oxygen-steam. ‘The atomic oxygen spectrum weakened with increasing 
dilution by steam, but a few lines remained even with pure steam. The hydroxyl- 
band and Balmer spectra were prominent throughout. 

Mercuric oxide films were formed, as with dry oxygen, even in mixtures con- 
taining only 5 per cent of oxygen. The ash-coated probes incandesced in all 
mixtures containing up to 95 per cent oxygen, but were brightest in steam-rich 
mixtures. 

(vil) Neon-steam. 'The hydroxyl-band and Balmer spectra only gave way 
to the neon spectrum as the neon concentration approached 100 per cent. The 
behaviour of the ash-coated and wire probes was similar to that observed in the 
case of the corresponding oxygen-steam mixtures, in that the incandescence 
effect increased with increasing steam concentration. 

The following conclusions have been drawn: 

Dry molecular hydrogen is not dissociated but only excited in the ring 
discharge. ‘The generally accepted view, according to which atomic hydrogen 
is formed and the absence of the Balmer spectrum is due to recombination of the 
atoms on the walls cf the discharge vessel, is therefore untenable. 

Atomic hydrogen issues in considerable quantity from the discharge in steam, 
hydrogen—steam or hydrogen—oxygen mixtures. The incandescence of ash-coated 
probes, and the spontaneous glowing of copper, tungsten, or platinum wires in 
the gases drawn from the discharge testify to a high concentration of atomic 
hydrogen, whilst the reduction of mercuric oxide films serves as a more sensitive 
test. It is unlikely that moisture adsorbed on the walls of the discharge vessel 
materially affects the yield of atomic hydrogen ; for, if it did, it should also cloak 
the catalytic activity of a thin copper film, which it does not do. 

The following reactions take place in hydrogen-steam and the equivalent 
hydrogen—oxygen mixtures. Steam is dissociated into hydroxyl and atomic 
hydrogen. The reverse reaction cannot proceed except through three-body 
collisions, and is therefore practically excluded. The same is true of the union 
of hydrogen atoms to molecules; consequently atomic hydrogen has a long 
survival life. On the other hand, the life of hydroxyl is short and is terminated 
by reacting with molecular hydrogen to form steam and atomic hydrogen, and with 
hydroxyl to form steam and atomic oxygen. Owing to the ease with which 
oxygen atoms form negative ions, they revert quickly to the molecular state. 
Finally, hydrogen reacts with oxygen to form steam and atomic oxygen. Some 
hydrogen peroxide may also be formed, but mainly at surfaces. 

Dry oxygen, unlike dry hydrogen, is dissociated directly into atomic oxygen, 
which has a short life and quickly reacts within the discharge to form molecular 
oxygen and ozone. With stagnant mercury vapour, ozone forms films of yellow 
mercuric oxide. All the reactions occurring in oxygen and in hydrogen-steam 
mixtures also take place in oxygen—-steam and the equivalent oxygen—hydrogen 
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mixtures ; thus ozone and atomic hydrogen appear together in the effluent gases 
and interact, if at all, very slowly. ‘The phenomena of the formation of mercuric 
oxide and the heating of ash-coated probes and wires in oxygen—steam mixtures 
are due to ozone and atomic hydrogen respectively, and are not evidence, as 
hitherto supposed, of a long life of oxygen in the atomic state. 

The products leaving the ring discharge fed with a mixture of steam (or the 
equivalent hydrogen—oxygen mixture) and sufficient mercury vapour consist 
almost whelly of steam ; only traces of atomic hydrogen are present. In this case 
the steam is no longer directly dissociated into hydroxyl and atomic hydrogen, but 
reacts with mercury to form hydroxyl and mercurous hydride. Other reactions 
occurring are the oxidation of mercurous hydride to form mercuric oxide and 
hydroxyl, the reduction of mercuric oxide by atomic (or molecular) hydrogen to 
mercury and hydroxyl (or steam), the reduction of mercury by molecular hydrogen 
to mercurous hydride and atomic hydrogen, and the formation of steam and atomic 
oxygen from hydroxyl. 

In neon-steam and the equivalent neon—hydrogen-oxygen mixtures, the 
reactions are similar to those occurring in steam. 
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The Significance of the Observation of Intense 
Radio-Frequency Emission from the Sun 
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ABSTRACT. A summarizing account is given of possible mechanisms which have been 
proposed to account for the emission of intense radio waves from the sun. Whilst most 
authors agree that the radiation from the undisturbed sun is due to an electron temperature 
of a million degrees in the solar corona, the greater intensity emitted from sunspots has led 
‘some authors to propose non-thermal mechanisms, based on the coherent oscillation of a 
large number of electrons in the solar corona above the sunspot. The mechanisms proposed 
for the maintenance of coherent oscillations are examined in detail in this paper and it is 
concluded that whilst they give adequate explanations for the oscillation observed in 
discharge tubes, no satisfactory mechanism has yet been suggested for the maintenance of 
electron oscillations in the solar corona. 

It is the1efore concluded that the observations can only be accounted for by the occur- 
rence, in the solar corona near sunspots, of electron temperatures of up to 10!° deg. k. 

The problem of maintaining an electron temperature of 10® deg. kK. in the normal corona, 
and 10% deg. k. in the corona near a large sunspot, has been considered in an earlier paper ; 
it was there shown that by the action of the general magnetic field of the sun and the non- 
uniform rotation of the solar surface, a potential difference is developed between the poles 
and the equator which is capable of maintaining an adequate electron temperature in both 


regions. 


§1. INTRODUCTION 


ECENT experimental work has shown that the sun emits radio waves of 
great intensity. For wavelengths between 1 metre and 5 metres, the 
intensity corresponds to black-body emission at temperatures of about 

10® deg. k. for the undisturbed sun and up to 10?° deg. K. for regions near sunspots. 
‘There is general agreement that the radiation from the undisturbed sun is due to 
a genuine electron temperature of about a million degrees in the solar corona. 

Some authors have suggested that the more intense radiation from regions 
near sunspots is produced by the coherent oscillation of a large number of electrons 
by some mechanism similar to that which produces plasma oscillations. If 
this were the true explanation it would not be necessary to assume that, near 
sunspots, the electrons had a mean energy corresponding to a temperature of 
10!° deg. K. 

The main object of this paper is to examine in detail the mechanisms which 
have been proposed for the coherent oscillation of electrons in the solar corona. 
It is shown that although it is a simple matter to demonstrate the occurrence in 
the solar corona of anatural “ plasma’”’ frequency of the electron gas, itis much more 
difficult to suggest any mechanism which would maintain the plasma oscillations, 
or any mechanism by which radio waves could be emitted from the plasma even 
if it did oscillate. All the mechanisms proposed to explain the oscillations in 
discharge tubes require sharp gradients of electron density and electron beams 
having a highly ordered motion; it is difficult to see how either of these conditions 
could occur'in the solar corona. The conclusion is therefore reached that no 
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satisfactory mechanism has yet been suggested. It is therefore concluded 
that the radiation emitted from near sunspots must be produced by the non- 
coherent acceleration of electrons, so that its observation implies the occurrence 
of genuine electron temperatures of 10!° deg. k. in these regions. A mechanism 
for the maintenance of such a great electron temperature in the neighbourhood 
of sunspots has been considered in an earlier paper. It was there suggested that 
by the action of the general magnetic field of the sun and the non-uniform rotation 
of the solar surface, a large potential difference is developed between the poles. 
and the equator, which is capable of maintaining a mean electron energy in the 
undisturbed corona corresponding to a temperature of at least 10°deg.k. The 
distortion cf the general magnetic field by the magnetic field of a sunspot allows. _ 
the maintenance of still greater electric fields in the corona near the sunspot, 
and mean electron energies corresponding to a temperature of 10'°deg.k. are 
likely to occur near large sunspots. 


§2. THE EXPERIMENTAL OBSERVATIONS 


The earliest deductions that the sun occasionally emits radio waves whose 
intensity greatly exceeds that expected from a source at 6,000°k. were made by 
Appleton (1945) and Hey (1946) from observations on frequencies of about 
70 Mc/s. Since that time a large amount of experimental work has been carried 
out on frequencies between 60 Mc/s. and 3,000 Mc/s. It has been found that 
for frequencies between 60 Mc/s. and 200 Mc/s. the intensity from the undisturbed 
sun corresponds to what would be emitted if the whole solar disc were at a tem- 
perature of about a million degrees (Pawsey 1946, Ryle and Vonberg 1948). 

At times of sunspot activity much greater intensities have been observed 
(Appleton 1945, Hey 1946, Pawsey, Payne-Scott and McCready 1946, Ryle and 
Vonberg 1946); the radiation at frequencies of about 200 Mc/s. then corresponds. 
to what would be emitted if the whole solar disc were at a temperature of about 
108 deg. K. 

Experiments have also been made to determine the position and size of the 
source of the intense radiation (Ryle and Vonberg 1946, McCready, Pawsey and 
Payne-Scott 1947). ‘The results obtained have shown that the radiation is emitted 
from a comparatively small region in the neighbourhood of the sunspot. The 
observed intensities of radiation therefore correspond to a still greater effective 
temperature, and it must be concluded that regions near sunspots are capable 
of emitting radiation whose intensity corresponds to the emission from a black- 
body source at a temperature exceeding 101° deg. k. 

The intense radiation from sunspots has been found to be circularly polarized, 
but the radiation from the undisturbed sun is apparently unpolarized (Martyn 
1946a, Appleton and Hey 1946, Ryle and Vonberg 1946, 1948). 


§3. EQUILIBRIUM MECHANISMS LEADING TO THE DEDUCTION 
OF VERY HIGH TEMPERATURES 

Some theories of radio-frequency emission from the sun have been based upom 
“equilibrium mechanisms”’, in which the intensity of the radiation tends to the 
value associated with the mean electron energy. The radiation from an infinite 
heated ionized gas of low density would correspond to that from a black-body, 
although in practical cases the total ‘‘ optical depth”? may be insufficient for the 
full intensity to be established. It is clear that such theories cannot account 
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for radiation of intensity greater than that emitted by a black-body having the 
same temperature as that of the electrons; an explanation of solar radiation by 
these theories therefore requires the existence of very high electron temperatures. 


(i) Undisturbed Sun 


In a previous communication (Ryle 1948) a mechanism for the generation of 
radio-frequency radiation in the solar atmosphere was developed. By an applica- 
tion of Appleton’s magneto-ionic theory to the propagation of metre waves in 
the solar atmosphere, it was shown that radiation could only escape to the earth 
from a source situated at or above the level at which the refractive index vanished 
(about 5 x 10®°cm. above the photosphere for a wave of frequency 200 Mc/s.). 
It was further shown that the absorption of an incoming wave was negligible 

except near the level at which the refractive index vanished, but that the total 
absorption in this region was considerable. It was therefore concluded that each 
region in the solar atmosphere would radiate at a particular frequency with an 
intensity which might approach that corresponding to a black-body radiator 
having the same temperature as the electrons in the region. It was therefore 
‘suggested that the observed intensity of the radiation from the undisturbed sun 
was the result of an electron temperature of at least a million degrees in the corona. 

It was suggested that the maintenance of such an electron temperature could 
be accounted for by the existence of a potential difference between the poles and 
equator of the sun, the potential difference being developed by the action of the 
general solar magnetic field on the non-uniform rotation of the photosphere. 

Martyn (1946b, 1948) has also developed a “‘thermal”’ theory of the radiation 
from the undisturbed sun, and whilst not suggesting a mechanism for the main- 
tenance of the electron temperature, he considers that the observations can be 
accounted for by assuming a temperature of a million degrees in the corona. 
Waldmeier (1948) and Shklovsky (1947) have also developed theories involving a 
comparable coronal temperature. Moreover, the observation of lines of very 
highly ionized atoms in the coronal spectrum leaves little doubt that the radiation 
from the undisturbed sun is due to purely “‘ thermal” processes in the corona. 


(ii) Sunspot Radiation 


In the previous paper (Ryle 1948) the theory was extended to the case of sunspot 
radiation. It was shown that in the solar corona above a sunspot appreciable 
absorption could occur at two levels (about 10cm. and 10!°cm. above the photo- 
sphere for a frequency of 200 Mc/s. and a sunspot magnetic field of 3,000 gauss) ; 
each region would absorb one of the circularly polarized components. It was 
therefore suggested that circularly polarized radiation would be generated in 
each of these regions; the intensities of the two waves would approach the values 
associated with the electron temperatures in the two regions. 

It was further shown that the distortion of the general solar magnetic field 
by the magnetic field of the sunspot could give rise to a very much greater electric 
field in the corona above a sunspot than in the undisturbed corona. It was 
shown that if the non-uniform rotation of the sun produced a potential difference 
between poles and equator which was capable of maintaining a normal electron 
temperature of 10®deg.x., then the distortion of the magnetic field near a large 
sunspot might be expected to lead to temperatures of the order of 10° deg. x. 
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§4. MECHANISMS INVOLVING THE COHERENT OSCILLATION OF 
ELECTRONS 


(i) Non-thermal Theories of Sunspot Radiation 


The extremely great electron temperatures which are required to account 
for the observed intensity of the radiation from sunspots have led other authors. 
to propose non-thermal mechanisms. In order to produce an intense radiated 
field by electrons of comparatively low energy, it seems necessary to develop a 
mechanism for maintaining the coherent oscillation of a large number of electrons. 

Martyn (1947) and Shklovsky (1947) have both suggested that the coherent 
oscillation of electrons in the corona could be produced by a mechanism analogous 
to the plasma oscillations familiar in gas discharge tubes. Ifsucha process were to 
occur, the intensity of the radiation emitted could be very much greater than that 
associated with the thermal motion of the electrons in the corona. It is therefore 
important to examine whether oscillations of this kind could occur under the 
conditions likely to exist in the solar corona. 


(ii) Existing Theories of the Mechanism of Plasma Oscillations 


Before attempting to consider the possibility that plasma oscillations could 
account for the intense emission from the solar corona it is important to understand 
the precise mechanisms which are envisaged in such theories. Unfortunately, 
the literature of plasma oscillations is notable for the small amount of detailed 
theoretical study as compared with the quantity of careful experimental results 
which have been obtained in the study of gas discharges. 

In the original theory of Tonks and Langmuir (1929) it was shown that the 
electrons in an ionized gas had a natural frequency of oscillation given by 
(Ne?/zm)?, where N is the electron density and e and m the charge and mass 
of the electron. (They also showed the existence of ionic oscillations having 
much lower frequencies, but it will be seen later that such oscillations cannot 
account for the emission of radio-frequency radiation from the sun.) It should 
be noted that although this theory provided natural periods of oscillation of the 
same order as those observed in experiments on gas discharges, no definite 
mechanism was proposed for the excitation of such oscillations. Without some 
additional mechanism, a plasma is not capable of generating oscillations, and 
has merely the characteristics of a tuned circuit. If some feedback mechanism 
or “‘negative resistance”’ is introduced, however, a plasma will then become 
capable of generating oscillations. 

A further difficulty arises in accounting for the detection of the oscillations. 
This difficulty is due to the fact that the group velocity in the plasma is zero 
(Tonks and Langmuir 1929, Tonks 1929); it is therefore impossible for energy 
to be propagated by an electromagnetic field from the oscillating region to the 
measuring probe, unless the gradient of electron density or magnetic field is. 
sufficiently great to allow appreciable transfer by way of the evanescent wave. 

The original theory has since been extended to account both for the excitation 
of plasma oscillations and for the escape of energy from the oscillating region. 


(a) The excitation of plasma oscillations. 


Later work by Merrill and Webb (1939) suggested that the oscillation of a 
thin layer near the cathode was responsible for velocity modulation of the primary 
electron stream from the cathode, a conclusion which has been confirmed by 
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Armstrong (1947). Recently Neill (1949) has shown that it is possible to account 
for the experimental observations by a feedback mechanism analogous to that 
by which a diode can maintain oscillations in virtue of transit-time effects. He 
supposes that electrons emitted by the cathode are subject to an accelerating 
field which depends on the space-charge in the plasma-sheath; if this space-charge 
is varying, the electron stream will be velocity-modulated, and if the distance 
to the plasma is suitable it will produce a bunching of the electrons which in 
turn can affect the space-charge. In this’way part of the energy from the cathode 
fall of potential is converted into energy of oscillation of the plasma. Neill 
has shown that the theory provides good quantitative agreement with many 
of the related observed phenomena. 

A more elaborate mechanism for the generation of oscillations in vacuum 
tubes has recently been suggested by Haeff (1949), which depends upon the 
interaction of two coherent beams of electrons.* Haeff has shown that under 
certain conditions a region containing two interpenetrating electron beams is 
capable of amplifying an initial disturbance having a frequency related to the 
plasma frequencies of the two beams. Whilst insufficient detail has been given 
of the precise mechanism, it appears that it is analogous to that of the travelling- 
wave tube, in that the disturbance is guided along the moving electron beams, 
and that the wave frequency to an observer moving with the electron beam is 
related to the plasma frequency. 

This brief survey shows that all the detailed mechanisms so far proposed for 
the maintenance of oscillations in a gas discharge tube depend on the velocity 
modulation of an incident electron beam. It will be necessary to examine 
(in § 4 (iii)) whether such a mechanism could provide excitation for the coherent 
oscillation of electrons in the solar corona. 


(b) The escape of energy from the oscillating region. 


Tonks and Langmuir (1929) considered that even in the small-scale conditions 
of the discharge tube, where the entire plasma occupies a thickness considerably 
less than a vacuum wavelength, it might be necessary to account for the detection 
of the oscillations by virtue of the escape of accelerated electrons from the plasma 
region. Merrill and Webb (1939) concluded that the alternating voltage observed 
on the measuring probe in their experiments was due to the periodic fluctuations 
of the space-charge in the incident electron stream. Later experiments by 
Armstrong (1947) and Neill (1949) have shown that the oscillating plasma will 
cause velocity modulation of the incident electron stream; they have also confirmed 
the earlier conclusion that the voltage observed on the measuring probe is caused 
by the bunching cf such an incident electron stream. 

Velocity modulation of the incident stream of electrons could also account 
for the escape of energy from a region which was oscillating by the mechanism 
suggested by Haeff (1949). The possibility that electromagnetic radiation could 
escape from such an oscillating region cannct be excluded until a more precise 
physical picture of the mechanism has been given; since, however, the emission 
of an electromagnetic wave relative to a coordinate system moving with the stream 
of electrons corresponds closely to the escape of radiation from a normal plasma, 
it seems probable that the escape of energy by electromagnetic radiation presents 
considerable difficulty. 


* A mechanism somewhat analogous to that proposed by Pierce (1948) to account for low- 
frequency oscillations in an electron stream passing through a stationary ion cloud. 
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In the absence of further information it therefore appears that energy is only 
likely to escape from an oscillating plasma by virtue of the velocity modulation 
imposed on an electron stream which has passed through the region. 

It is clear that in the large-scale conditions of the solar corona, the escape 
of energy from an oscillating plasma presents even more serious difficulties 
than for the case of the discharge tube; the application of the theories outlined 
above to the solar corona will be discussed in the next section. 


(iti) The Possibility of Electronic Plasma Oscillations in the Solar Corona 


It is first important to note that the existing experimental and theoretical 
work on plasma oscillations in gas discharge tubes relates to conditions in which 
the dimensions of the plasma are small compared with the vacuum wavelength 
of the oscillation. In the solar corona it is probable that the electron density 
and magnetic field will be virtually constant over a volume of many cubic 
wavelengths. 


(a) The excitation of plasma oscillations in the corona. 


In a gas discharge tube containing a cathode at a low temperature, the initial 
velocities of the electrons are small compared with the velocity produced by the 
applied electric field. Under the influence of a steady applied field, the stream 
of electrons will therefore travel a distance of the order of one mean free path 
before the motion becomes disordered. If now an alternating field is super- 
imposed on the steady field near the cathode, the electrons will be velocity- 
modulated and, if the mean free path is sufficient, will become bunched at some 
distance from the cathode. Neill (1949) has shown that the variation of the space- 
charge resulting from this bunching can then provide excitation to the oscillating 
plasma. 

If now the temperature of the cathode is increased, the initial distribution 
of the velocities of the electrons may become comparable with the velocity produced 
by the applied field; although the superimposition of an alternating fiela will 
still result in the acceleration of each electron, the spread in the initial velocities 
of the electrons will prevent effective bunching. An alternating density of space- 
charge will therefore no longer be produced, and the oscillations of the plasma 
cannot be maintained. 

In the solar corona the initial random velocities of the electrons in any particular 
region will be of the same order as the velocities gained during a mean free path, 
and the conditions therefore correspond to those in a discharge tube having a 
cathode at a very high temperature. The superimposition of an alternating 
field will not produce appreciable bunching in such circumstances ; the mechanism 
suggested by Neill (1949) cannot therefore account for the maintenance of plasma 
oscillations in the solar corona. 

Although Haeff (1949) has given no detailed description of the conditions 
which are necessary for the generation of oscillations by the interaction of two 
electron beams, it seems likely that here also the distribution of velocities in the 
primary electron beams must initially be small compared with the ordered velocity. 

It appears that any mechanism which involves velocity modulation of an 
incident electron stream is equally unlikely to be effective in maintaining the 
coherent oscillation of electrons in the conditions of the solar corona. 
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(b) The escape of the energy. 


Even if itis assumed that the coherent oscillation of electrons could be maintained 
in the solar corona, it is doubtful whether any energy so generated could escape 
‘to the earth. It is clearly more difficult to account for leakage of the evanescent 
wave in a region where the gradient of electron density and magnetic field will 
be very much smaller than those in a discharge tube. It is therefore necessary 
‘to examine the possibility that energy might be transferred from the oscillating 
region by the passage of electrons through the region. A stream of electrons 
undergoing velocity modulation in this way might pass into an overlying region 
‘from which radiation could escape; such a stream is unlikely to be an efficient 
radiator, but some of the energy imposed by the oscillating plasma-might be able 
to escape. It is, however, necessary that the electron stream shall still be effec- 
‘tively bunched when it emerges into the region from which the radiation can be 
‘propagated. Unless there are very considerable gradients of electron density 
-or magnetic field, the electron stream must travel a distance corresponding to a 
“very large number of wavelengths. The accurate collimation of a velocity- 
modulated electron stream which has traversed such a distance can only occur 
-if the random velocity of the electrons is negligible compared with their ordered 
motion. It has already been shown that such a condition is unlikely in the solar 
‘corona. 

It therefore appears that although the electrons in the solar corona may have 
-a Suitable natural period of oscillation, there are serious difficulties in accounting 
‘for the intense radiation at radio frequencies by mechanisms involving plasma 
oscillations. Satisfactory theories have been developed to account for both the 
-driving mechanism and the mechanism of energy transfer for plasma oscillations 
‘in discharge tubes; these theories are not applicable in the large-scale conditions 
-of the solar corona unless very large gradients of electron density or magnetic 
‘field are assumed. ‘The maintenance of such gradients in a region where the mean 
‘free path for an electron is many kilometres does not appear possible. 


(iv) The Possibility of Ionic Plasma Oscillations in the Solar Corona 


So far attention has only been given to the coherent oscillation of electrons. 
If positive ions are also present, they will be capable of oscillating with a natural 
‘frequency of (Ne?/7M)?, where N and M are the number density and mass of the 
ion and e is the electronic charge (Tonks and Langmuir 1929). In a region 
containing an equal number of ions and electrons, the natural frequency of the 
ionic oscillations will be very much less than that of the electronic oscillations ; 
‘the refractive index for a wave having the frequency of the ionic oscillations is 
therefore imaginary, and ionic oscillations cannot generate radio-frequency 
radiation which can escape through the overlying layers. 

Tonks and Langmuir (1929) have, however, also shown that oscillations of 
-very low frequency can be generated in an ionized gas, and that these oscillations 
may be propagated in a manner analogous to the propagation of sound waves. 
It is possible that such oscillations might occur in the solar corona; the passage 
of such a wave through a region which was emitting radio-frequency radiation 
would produce variations in the mean energy or density of the particles, and might 
-cause fluctuations in the intensity of the emission. ‘Thus, whilst ionic plasma 
oscillations cannot account for the great intensity of the radio-frequency radiation 
-emitted by the corona, it is possible that they might produce low-frequency 
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fluctuations in the intensity of the emitted radiation. Such waves might be 
responsible for the “bursts” of radiation associated with the emission from 
sunspot regions. 

§5. CONCLUSIONS 

A detailed consideration of theories of plasma oscillations has not revealed 
a satisfactory explanation for the maintenance of coherent electronic oscillations 
in the solar corona. It is therefore concluded that the observed radio-frequency 
radiation arises from purely random motion of the electrons in the corona. It 
has previously been shown (Ryle 1948) that each region of the corona is absorbing 
to a certain radio frequency, and thus each region will radiate at a particular fre- 
quency with an intensity depending upon the electron temperature and the total 
absorption coefficient. If the absorption is sufficient, the intensity will equal 
that of black-body radiation corresponding to the electron temperature in the 
region. 

If these conclusions are accepted, the observations show that the electron 
temperature in the normal corona is at least a million degrees (a figure assumed 
in many theories of the radiation from the undisturbed sun), whilst in the corona 
above large sunspots the electron temperature may reach 10!°deg.x. The 
existence of such electron temperatures in the solar corona is clearly of great 
significance; observations over a range of radio frequencies may allow a detailed 
study to be made of the distribution of electron temperature throughout the 
corona. 

A mechanism for the maintenance of a large electron temperature in the solar 
corona has already been suggested (Ryle 1948). It was shown that the non- 
uniform rotation of the solar surface, together with the general magnetic field of 
the sun, will produce a potential difference of about 5 x 10° volts between poles 
and equator. A potential difference of this order is sufficient to maintain not 
only a normal electron temperature in the corona of about a million degrees, but 
also an increased electron temperature in the corona above sunspots. For a 
large sunspot the potential gradient is sufficient to maintain in the regions 
responsible for the emission of radio-frequency radiation a mean electron energy 
corresponding to a temperature of 10! deg. k. 

The application of similar mechanisms to account for the discrete sources of 
radio waves in the galaxy, and the possible relation between these sources and the 
origin of cosmic rays, will be considered in another paper (Ryle 1949). 
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ABSTRACT. Several authors have suggested the possibility that cosmic rays are due to 
the acceleration of charged particles in the atmospheres of certain stars. Recent obser- 
vations at radio frequencies have provided some evidence in support of this hypothesis 
which is reviewed in this paper. 

The observations have shown that at least part of the radio-frequency radiation from the 
galaxy is emitted by discrete sources of small angular diameter. Measurements of the 
intensity of the radiation from these sources, and of the time variations of the intensity, 
indicate that the sources emit radio waves as if they were at a temperature of more than 
1014 deg. k. It is concluded that this intense radiation cannot be caused by the coherent 
oscillation of a large number of electrons, but must be due to a genuine electron temperature 
of about 10 deg. kK. (corresponding to a mean electron energy of 101° electron volts). The 
existence of this electron temperature indicates that in certain stellar bodies there are 
mechanisms capable of accelerating particles to cosmic ray energies. 

A previous theory of the emission of radio waves from the sun is extended to other 
stellar bodies. It is concluded that a star having a surface magnetic field strength and a 
peripheral velocity somewhat greater than those which have so far been observed could emit 
radio waves having the observed intensity. A star of this type could also accelerate charged 
particles to cosmic ray energies. 

The small visual brightness of the sources of radio waves suggests that these bodies may 
be characterized by small visual opacity or small photospheric temperature. 

It is suggested that the difference between the angular distribution of cosmic rays and of 
the radio waves from the galaxy is due to the deflection of cosmic ray particles by magnetic 
fields. In addition to previous theories of a general magnetic field in the galaxy, it is shown 
that the scattering produced by the magnetic fields of the distributed sources themselves 
should be sufficient to produce an isotropic distribution of cosmic rays at the earth even 
though the sources may show a marked concentration near the galactic equator. 


§1. INTRODUCTION 

EVERAL authors have proposed mechanisms for the generation of cosmic 

rays by the acceleration of charged particles in the atmospheres of certain 

stars. Swann (1933) demonstrated that the growth of the magnetic field 

of a sunspot could give rise to potentials of the order of 10° volts in.a region where 
little energy would be lost by collisions. He suggested that similar disturbances 
in other stars might give rise to particles having an energy of more than 101° 
electron volts. Alfvén (1937) showed that even in the absence of sunspots a 
comparable potential could be developed by the non-uniform rotation of the solar 
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surface. Davis (1947) deduced that the star BD-18°-3789, which recent measure- 
ments (Babcock 1947b) have shown to have a magnetic field strength at the 
poles of about 5,500 gauss, would produce, relative to a stationary observer, a 
potential difference of about 10!*volts. He showed, however, that owing to 
the formation of a space-charge this potential difference was not available for the 
acceleration of charged particles. Babcock (1948) has considered this process 
in relation to magnetically variable stars, and he has shown that the periodic 
reversal of the magnetic field would allow the escape of particles of high energy 
during a certain fraction of the magnetic cycle. 

Recent measurements (Hey, Parsons and Phillips 1946, Bolton and Stanley 
1948, Ryle and Smith 1948) on the emission of radio-frequency radiation from the 
galaxy have provided independent evidence for the existence of such sources; 
it is the purpose of this paper to examine this evidence and to make certain deduc- 
tions concerning the nature of the sources. 

The observations have shown that at least part of the radio-frequency radiation 
from the galaxy is emitted by discrete sources, whose angular diameters are smaller 
than the resolving power of the apparatus (about eight minutes of arc). It has 
also been found that the intensity from some of the sources fluctuates rapidly 
and that sometimes ‘‘bursts’’ of radiation are observed, in which the intensity 
is doubled for a period of about 20 seconds. Providing that these fluctuations 
represent genuine variations of the emission from the source, it is possible to fix 
an upper limit to the dimensions of the source. Extended measurements of 
the positions of the two intense sources have suggested that they are not part 
of the solar system; if it is assumed that they are situated at the distance of the 
nearest fixed star, the deduced diameter and the observed intensity would corre- 
spond to the emission from a black-body source at a temperature of more than 
10 der x: 

It is possible that some special radiating process might be responsible for the 
emission of radiation having an intensity greater than that corresponding to the 
electron temperature; in the absence of such a process, the radiation can only be 
accounted for by assuming that the electrons have a random energy corresponding 
to the thermal energy at a temperature of 101deg.x. (the mean energy of an 
electron is then about 10?° electron volts). 

The maintenance of this electron energy implies the existence of mechanisms 
capable of accelerating particles to cosmic ray energies. It is therefore suggested 
that the sources of intense radio-frequency radiation may also be responsible for 
the emission of cosmic rays. 

The assumptions made in the development of this theory are analysed in 
detail in the following sections. A mechanism is suggested by which intense 
radio waves (and cosmic ray particles) could be emitted from the atmospheres 
of certain stars having a large rotational velocity and surface magnetic field. 
It is further shown that the scattering effect caused by the magnetic fields of such 
stars would be expected to lead to an isotropic distribution of cosmic rays at the 
earth, independent of the distribution of the sources in the galaxy. 


§2. THE EXPERIMENTAL EVIDENCE 


Hey, Parsons and Phillips (1946) first observed that the intensity of the radio- 
frequency radiation from the constellation of Cygnus showed fluctuations in a 
period of a few minutes; they therefore concluded that at least part of the radiation 
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from the galaxy was due to the emission from discrete sources. More refined. 
experiments by Bolton and Stanley (1948) and Bolton (1948) enabled the source 
to be located more precisely, and showed that its angular diameter was less than. 
the resolving power of the apparatus (about eight minutes of arc). Similar 
observations by Ryle and Smith (1948) showed a second, more intense, source 
in the constellation of Cassiopeia, and extended observations showed that whilst. 
the intensity from the source in Cygnus frequently showed sudden increases, 
that from the source in Cassiopeia was notably constant. They also found that 
the radiation from both sources appeared to be randomly polarized. 

The intensity of the 80 Mc/s. radiation from the sources in Cygnus and 
Cassiopeia has been found to be of the order of 2 x 10-?2w.m-?(c/s.)-1. Both 
Bolton (1948) and Ryle and Smith (1948) have detected other discrete sources. 
of smaller intensity, and some of these have also shown fluctuations of intensity.. 
None of the sources appears to coincide with an outstanding visual star. 


§3. DETERMINATION OF THE EQUIVALENT TEMPERATURE 
OFA SOURCE 

The intensity of the radiation emitted by a distant source may be compared 
with that from a black-body source which subtends an equal solid angle. If 
the intensity corresponds to that produced by the black-body at a temperature 7, 
the source may be said to have an “‘ equivalent temperature” JT. The significance 
of this temperature will be discussed in §4, and it will be used now merely as a 
convenient unit with which to express the emission from the source. 

In order to determine the equivalent temperature of a galactic source of radio- 
frequency radiation, it is necessary to observe both the intensity of the radiation 
at the earth and the solid angle subtended by the source. Unfortunately it has 
not been possible to achieve a sufficiently large resolving power to measure the 
solid angle directly. Observations of the source in Cygnus have, however, shown 
sudden fluctuations in which the intensity is doubled in a period of 20 seconds. 
This variation might represent variations in emission, or might be caused by 
refraction in an intervening medium; it will first be shown that a refraction process 
is unlikely. 

The observation that the intensity from the source in Cygnus exhibits rapid. 
fluctuations, whilst that from the source in Cassiopeia remains noticeably constant, 
suggests that the fluctuations are not produced by general refraction in interstellar 
matter. Since the observations were made when the sources were nearly vertically 
overhead, it seems unlikely that refraction in the ionosphere could affect one source 
and not the other (the angles of incidence during observation are about 12° 
for the source in Cygnus and 6° for the source in Cassiopeia). Whilst the possi-- 
bility still remains that the intensity from the source in Cygnus is affected by 
refraction in some localized region of interstellar matter, it seems most probable- 
that the observed fluctuations represent a true variation of the emission. 

The total intensity of the radiation received from a source may be considered 
as being due to the addition of the radiation emitted by a large number of small. 
elements of the source. Ifthe radiation from each of these small elements exhibits 
sudden increases of intensity, the total intensity will only show such increases 
if the fluctuations from different elements arrive simultaneously. If large increases. 
are observed it is therefore necessary that the different elements of the source. 
shall be separated by a distance not greater than the distance travelled by an. 
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electromagnetic wave in the duration of the fluctuations. In the absence of any 
special orientation of the emitting surfaces relative to the earth it is therefore 
possible to deduce a maximum dimension of the source by determining the duration 
of the sudden fluctuations. 

Observations of the source in Cygaus have provided records showing a doubling 
of the intensity in a time of 20 seconds. It is therefore deduced that the source 
of the variable component of the radiation occupies a volume whose linear dimen- 
sions do not greatly exceed the distance travelled by electromagnetic waves in 
_20 seconds (about 10!?cm.). 

There is, as yet, insufficient information on the distance of the sources. 
Parallax measurements at radio frequencies are unlikely to be of sufficient accuracy 
to determine the distances, although preliminary measurements by Ryle and 
Smith (unpublished) have shown that the sources in Cygnus and Cassiopeia 
are at distances greater than 2x 10!®cm. Improvements in the angular accuracy 

-of the measuring apparatus are unlikely to lead to a significant extension of parallax 
measurements by radio methods, although they may allow the sources to be 
identified with known stellar bodies, in which case parallax measurements by 
visual methods may be possible. The existing measurements indicate that it 
is unlikely that the sources are part of the solar system, and in the absence of 
further experimental evidence it will be assumed that they are not nearer than the 
‘nearest known fixed star (about 3 x 10!8cm.). 

With this assumption, and for the deduced maximum diameter of 10!¥cm., 
it is found that the source in Cygnus emits radiation at a frequency of 80 Mc/s. 
_as if it were a black-body at a temperature of 10" deg. K. 


§4. THE EVIDENCE FOR THE EXISTENCE OF HIGH-ENERGY 
PARTICLES IN THE SOURCE 


The occurrence in the source of a true temperature of 10'deg.K. would 
involve collisions by electrons having a mean energy of about 10! electron volts. 
It is necessary to examine the possibility that special mechanisms may be 
responsible for the emission of intense radio waves by electrons of comparatively 
low energy. 

The coherent oscillation of a large number of electrons might produce intense 
radio waves; both Martyn (1947) and Shklovsky (1947) have suggested that 
‘plasma oscillations in the solar corona could, in this way, account for the intense 
radiation associated with sunspots. ‘This mechanism has been examined in 
some detail in a previous communication (Ryle 1949) which will be referred to 
as Paper 2. It was there concluded that the mechanism could not account for 
the emission of intense radio waves from the sun. It is now necessary to apply the 
arguments then used to determine whether plasma oscillations could account 
for the intense radiation from galactic sources. 

When considering the possibility of plasma oscillations in the solar corona 
‘the range of electron density and magnetic field likely to occur in the region 
of the source were known from well-established observations. For the cease 
-of the galactic sources there is no information on the magnitude of the electron 
density or the magnetic field except that obtained from the radio-frequency 
observations. ; 

It is first important to note that the radiation will only be able to reach the 
-earth if the refractive index remains real at all points between the source and the 
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earth. In the absence of a magnetic field the refractive index for a wave having 
a frequency of 80 Mc/s. vanishes at an electron density of about 8 x 10’ cm-® 
In the presence of an intense magnetic field the propagation of one circularly 
polarized component becomes possible for an appreciably greater electron density. 
If the source were situated in a region of great electron density and an intense 
magnetic field permitted the escape of one component, the radiation observed 
at the earth would, however, appear circularly polarized, a result which is not 
confirmed by experiment (Ryle and Smith 1948). 

It is therefore concluded that the presence of an intense magnetic field is 
not essential to the emission of radio waves from the type of galactic source which 
has been described. (The variable component of the radiation from some 
galactic sources is suggestive of the radiation associated with sunspots, although 
preliminary experiments have suggested that this component is also randomly 
polarized.) 

It might be suggested that circularly polarized radiation emitted from the 
two hemispheres of a star would be of opposite polarity and approximately equal 
intensity for an observer situated near the equatorial plane of the star. The observer 
would therefore deduce that the radiation was randomly polarized. Detailed 
observations of the polarization from both the intense galactic sources have, 
however, shown that the polarized component appears to be less than 5% of the 
total intensity; it is unlikely that both sources would be orientated with such 
precision as to give this kind of signal at the earth if the radiation were comprised 
of two circularly polarized waves of opposite polarity. It is therefore concluded 
that the observed random polarization is genuine; both components can therefore 
escape from the region of the source. 

In the absence of an intense magnetic field in the region of generation, the 
radiation must be emitted at or above the level at which the electron density is 
about 8x10’cm~*. Unless the temperature is less than about 10,000°K., 
ionization will be complete and the total particle density will not greatly exceed 
BOSCO? cm-*. 

It is thus apparent that the observed radiation can only be generated in regions 
where the electron density and magnetic field are similar to those in the solar 
corona. In Paper 2 reasons were given for supposing that plasma oscillations 
could not occur under the conditions of the solar corona, and it was suggested 
that the observed intensity of radio waves was due to the random motion of electrons 
having a mean energy corresponding to a temperature of 10° deg. kK. in the normal 
corona and up to 101!°deg.k. near sunspots. ‘The same arguments may now 
be applied to the radiation from galactic sources; in the absence of coherent 
electronic oscillations the observed intensity can only be accounted for by an 
electron temperature of at least 1014 deg. kK. 

The emission of radio waves from an electron gas and the maintenance of a 
large electron temperature in the solar corona by an electric field have been 
considered in an earlier paper (Ryle 1948) ; this paper will be referred to as Paper 1. 
It seems probable that similar mechanisms occur in the sources of intense galactic 
radiation; if an electric field is responsible for the maintenance of an electron 
temperature of 101*deg. K., it is necessary that the energy gained by an electron 
between collisions shall be equal to the mean energy of thermal agitation of a gas 
at a temperature of 1014deg.x. ‘The electron density in the region of generation 
has been deduced from a determination of the refractive index at a frequency of 
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80 Mc/s.; the mean free path of an electron (for appreciable loss of energy) is. 
then about 10!°cm. It may therefore be seen that an electric field of 1 volt cm~? 
would be sufficient to maintain a mean electron energy of 10!° volts (corresponding 
to an electron temperature of 10! deg. K.). 

The existence of an electric field of this magnitude would be capable of 
accelerating charges of either sign to a very great energy and, as has been shown 
for the case of the solar corona (Paper 1), some of these particles would be able 
to escape from the stellar atmosphere with a large energy. 

It is therefore suggested that the observed galactic sources are also capable 
of emitting particles having an energy corresponding to a potential of at least 
10!° electron volts and that the total emission of particles from these bodies may 
be the origin of the primary cosmic rays. 

The application of a theory of Alfvén (1937) to account for the maintenance 
of a large electron temperature in the solar corona has already been considered 
in Paper 1; the extension of these mechanisms to account for greater electron 
temperatures in the envelopes of other stellar bodies will now be considered. 


§5. THE PRODUCTION OF HIGH-ENERGY PARTICLES 
IN A STELLAR ENVELOPE 

Alfvén (1937) has shown that the potential developed between the poles and 
equator of a star of radius Ry and polar field strength Hp, is of the order of 
H,R?(w —w') where w and w’ represent the angular velocities at the equator and 
the poles. For a star embedded in a region containing a large density of inter- 
stellar matter, the available potential may be significantly greater and is given by 
H, Rew. ‘The recent observation (Babcock 1947a) of stars ‘having magnetic 
field strengths and rotational velocities very much greater than those of the sun 
therefore suggests that in certain stars a very large potential difference may be 
developed. 

It was shown in Paper 1 that the intensity of radio waves emitted at any 
frequency by a stellar envelope depends on the mean electron energy in the region 
where the refractive index for that frequency vanishes. This energy depends 
on (a) the electric field and (6) the electron density, and the latter is determined 
by the frequency of the wave,* so that the intensity at a given frequency is deter- 
mined solely by the magnitude of the electric field. 

For a star in which the available potential is given by H)R>w, the electric field. 
in the envelope will be of the order Hj)Ryw. Observations by Babcock (1947 b) 
have shown that certain stars are likely to have, at the photosphere, a polar field 
strength of 5,000 gauss and a peripheral velocity (Row) of 2x 10’cm.sec-!. The 
product of these two quantities is about 10* times the corresponding figure for 
the sun; a star having these characteristics might therefore be expected to emit 
radio waves corresponding to a temperature about 10* times that observed for 
the sun. It is possible that stars having somewhat greater magnetic fields and 
rotational velocities may exist, and it is suggested that the mechanism described 
would be adequate to account for the maintenance of an electron temperature: 
as great as 101 deg. kK. in the envelopes of certain stars. 


§6. THE ANGULAR DISTRIBUTION OF THE SOURCES 
Observations by Reber (1944) and Hey, Phillips and Parsons (1946) have: 


shown that the intensity of radio waves is greatest from directions near the galactic: 


* This statement assumes that the emitted wave is unpolarized. 
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equator, whereas the distribution of cosmic rays at the earth is highly isotropic. 
If the hypothesis put forward in this paper is correct, it is necessary to show how 
the observed isotropic distribution of cosmic rays could be produced during their 
passage from remote sources to the earth; it will be suggested that the important 
controlling factor is the magnetic fields of those stellar bodies which themselves 
produce the particles. 

Babcock (1948) has shown that the magnetic field of certain stars could deflect 
a cosmic ray particle at a considerable distance; it is possible to calculate the 
“effective radius for large-angle scattering” for particles of different energies. 
If it is assumed that the stars responsible for scattering cosmic ray particles have, 
for particles of energy 10! electron volts, an “effective scattering radius” R, 
and are distributed through the galaxy with density A, then if p is the effective 
radius of the galaxy, the probability P that a particle is deflected through a large 
angle as it traverses the galaxy is given by 


PeakeAp: AF Sai Goel oot (1) 


If now it is assumed that certain stellar bodies are responsible for the production 
of the radio waves from the galaxy (and of the cosmic ray particles), and that these 
bodies also cause the deflection of the particles, it is possible to use the knowledge 
derived from the radio measurements to make an estimate of P. 

The average radius of these stellar bodies for the emission of radio waves will 
be of the same order as the ‘ photospheric”’ radius Ry; they will therefore subtend 
a total solid angle 7R2Ap at the earth. 

It has already been shown that the emission of radio waves at a frequency 
of 80 Mc/s. corresponds to an effective temperature of 10'*deg.k. at the source, 
so that if all the sources subtend a total solid angle 7R>Ap, the effective mean 
temperature measured at this frequency will be 7RjAp x 10'4deg.k. Measure- 
ments at 80 Mc/s. have shown that the mean temperature of the galaxy is about 
104 deg. k., and so it is deduced that 7R?Ap ~10-". 

Substitution from equation (1) now gives: 
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It would therefore be possible to deduce a value for P without a knowledge of 
A or p, providing (R/R,)? were known. 

From calculations based on Babcock’s theory (1948), it can be shown that for 
a cosmic ray particle having an energy of 10! electron volts, (R/ Ro)? =3 x 10-8HyRy 
(where electromagnetic C.G.S. units are used). Reasons have already been given 
(§5) for supposing that galactic radio waves are produced in bodies for which 
H,R,w is about 108 times the corresponding quantity for the sun. On the 
assumption that their angular velocity is the same as that of the sun, H)R)>=3 x 10° 
and (R/Ro)? ~ 1028. 

The probability that a cosmic ray particle of energy 101° electron volts will 
be deflected through a large angle in passing through the galaxy would therefore 
be 1,000 (from equation (2)). Even if the angular velocities of the stellar bodies 
were 1,000 times as great as that of the sun, the probability of a violent deflection 
would be unity and the distribution of cosmic rays at the earth would still be 
isotropic. 

It is therefore concluded that even in the absence of a general galactic magnetic 
field (Alfvén 1938, Babcock 1947b) considerable scattering would be produced 
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by the magnetic fields of the distributed sources themselves; the difference in 
angular distribution of cosmic rays and the radio waves from the galaxy is therefore 
not regarded as an argument against their generation in a common source. 


§7. CONCLUSIONS 


Observations of intense sources of radio-frequency radiation have suggested 
that certain stellar bodies are responsible for the emission, at metre wavelengths, 
of radiation whose intensity corresponds to the radiation from a black-body 
source at a temperature of 10/4deg.xK. This estimate is based on deductions 
made from the intensity and its time fluctuations and on the assumption that the 
source is not nearer than the nearest known fixed star. So far it has not been possible 
to identify any of the sources with known stars, but it seems unlikely that a 
body having the characteristics deduced from the radio observations could pass 
unnoticed if it were at a comparatively small distance. 

The possibility that intense radio waves could be emitted by the coherent 
motion of electrons of comparatively low energy has been examined, but the 
necessary conditions seem improbable in a stellar envelope. It is therefore 
concluded that the observations indicate the existence of regions in which the 
mean electron energy is of the order of 10'° electron volts. ‘These regions would, 
under certain conditions, be capable of emitting particles having energies of this 
order of magnitude, and it is suggested that they are responsible for the production 
of primary cosmic rays. 

The difficulty of identifying sources of such outstanding characteristics with 
those of visual observations may imply the existence of a type of star in which a 
high temperature corona is associated with a comparatively small visual tempera- 
ture or opacity. This condition might arise in rapidly rotating stars having a 
large magnetic field. 

The observations which have already been made on discrete sources of radio- 
frequency radiation suggest that it may be possible to explain the total radiation 
from the galaxy on the basis of a number of such sources. It therefore appears 
that a type of star which has so far not been classified in astronomy at visual 
wavelengths may be responsible for two major astrophysical phenomena—the 
origin of cosmic rays and of the intense radio waves from the galaxy. 

The difference in the angular distribution of the radio waves and cosmic rays 
is accounted for by the deflecting effects of magnetic fields in the galaxy. In 
addition to previous theories of interstellar magnetic fields it is shown that the 
magnetic fields of the sources themselves are likely to cause an isotropic distribution 
of cosmic rays at the earth. 
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Pair-Production in Different Elements by y-Rays 


BYahe hee ROY, 
Wheatstone Physics Laboratory, King’s College, University of London 


‘Communicated by F. T. Randall ; MS. received 5th October 1948, and in amended form 
6th April 1949 


ABSTRACT. ‘The production of electron pairs on bombardment of Pb, Pt, Sn, Ag and 
Cu by y-rays from 200 mg. of Ra and its equilibrium products has been investigated. 
1,200 photographs were taken for each element and the yield was 100 pairs in Pb, 85 in Pt, 
‘72 in Sn, 63 in Ag and 50 in Cu. The distribution of energy of the positrons in Pb and Pt 
follows the shape of the theoretical curve of Bethe and Heitler, but the position of the 
maximum of the energy distribution curves in Sn, Ag and Cu shows a discrepancy with that 
tto be expected theoretically. The theoretical difference between the average energies of the 
positrons and electrons in the elements studied, when plotted against Z, is consistently 
higher than the experimental, though both can be approximately represented by straight 
lines. The angular distribution of the positrons and electrons with respect to the incident 
‘y-rays agrees with the theory. 


$1. INTRODUCTION 


HE creation of a pair of electrons of opposite charge in the Coulomb field 
of the nucleus by y-rays of energy greater than 2m)c? is a well-known 
phenomenon, and the cross section of this process has been calculated by 
Bethe and Heitler (1934) and by Jaeger and Hulme (1936). The experimental 
cross section for pair-production in Pb has also been studied, in particular by 
Chadwick, Blackett and Occhialini (1934), and it was found to agree quite well 
with the theoretical cross section as given by Bethe and Heitler. ‘The present 
experiments were undertaken to investigate (i) the energy spectra oi the positrons 
of the pairs created in Pb, Pt, Sn, Ag and Cu, (ii) the difference between the 
average energy of the positrons and electrons as a function of the atomic number 
of the elements studied, (iii) the angle between the positron and electron and the 
angles of emission with respect to the direction of the incident y-rays. 


§2. EXPERIMENTAL ARRANGEMENTS 


_ The experiments were carried out in a rubber-piston type cloud-chamber 
16-5cm. in diameter and 4cm. in height; the standard source of illumination, 
viz. five 100-volt, 100-watt lamps, flashed across 200 volts, was used. A magnetic 
field of 400 gauss, supplied by a pair of Helmholtz coils, was used to determine the 
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energy of the tracks. ‘The measurements of the angles and curvatures of the tracks. 
were taken after re-projecting the images from the two cameras. ‘To reduce 
inaccuracy arising from multiple scattering, only those tracks whose radius of 
curvature lay between certain limits were taken into account. The source used 
was y-rays from 200 mg. of Ra in equilibrium with its products, and, although 
the y-rays from Ra have several lines above 1 Mev., in these experiments observa- 
tions were restricted by energy consideration to those y-rays the energy of which 
was 2:2mev. The foil of the element to be studied was placed in the cloud- 
chamber, against the glass ring; facing the side where the foil was situated was 
a lead block 45 cm. in length, 12 cm. in width, and 10cm. in height, through which 
a hcle 5mm. in diameter had been drilled. This lead was necessary to prevent 
continuous ionization in the chamber due to y-rays. At the other end of the block, 
away from the cloud-chamber, there was placed on one side of the hole a thin- 
walled iron container with the source inside it ; on the other side, an electromagnet. 
When the expansion took place, the container was drawn towards the magnet; 
hence the source was in line with the hole and the y-rays could pass through and. 
bombard the element. 


§3. RESULTS AND DISCUSSION 


The results are given in Table 1. 


Table 1 
Element studied Pb RE Sn Ag Cu 
Thickness of foil () 15 15 20 20 25 
No. of photographs taken 1200 1200 1200 1200 1200 
No. of pairs found 100 85 72 63 50 


Taking into account the density and thickness of the material, the yield of pairs 
bears out, approximately, the theoretical predictions. 


(i) The Energy Spectrum of the Positrons in Pb, Pt, Sn, Ag and Cu 


As mentioned above, the energy of the incident y-rays was 2:2 Mev., and hence: 
the energy available for the pair is 1:-2mMev. From the experimental results it 
was found that the energy spectra for Pb and Pt, and for Ag and Sn, were essentially 
the same; therefore, only the energy distribution curves for Pb and Sn (Figures 
1 (a) and 1(6)), and that for Cu (Figure 1(c)), are given. They are plotted against 
the theoretical curve of Bethe and Heitler. 

The abscissa in each figure represents the total energy ot the positrons, that is,. 
the sum of the kinetic and rest energies ; the ordinate gives the ratio of the number 
of positrons at a particular energy to that of the total number of positrons observed 
from a given element. Examination of these curves shows that in the case of 
Pb and Pt the maximum lies at about 1-38 Mev., in the case of Sn and Agat 1:25 Mev. 
and in the case of Cu at 1-15 mev. It will be seen that for 2:2 Mev. incident y-rays. 
the energy distribution curve for Pb and Pt agrees in shape with the theoretical. 
curve. However, whereas the energy at which the maximum of the theoretical 
curve occurs remains the same for all elements, the experimental results show a 
variation in the maximum from higher to lower values according as the element 
studied is heavier or lighter. In addition, it may be noted that the curve of Pb 
and Pt falls very rapidly above 1-4 Mev. and there are more particles at the higher 
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‘than at the lower energies. The curve of Cu is almost symmetrical and the energy 

of the y-quantum after the creation of a pair is more or less equally distributed 
between the positron and electron. 

In order to calculate the energy of the positrons the curvatures of the tracks 
have to be measured. This measurement involves an inaccuracy at low energies, 
because of multiple scattering and the difficulty in ascertaining the small radii 
of curvature. At high energies inaccuracies are also liable to occur because of 
the large radius of curvature. The maximum error here was not more than 8°% 
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Figure 1 (a). Figure 1 (0). 


Figure 1.(c). 


in the extreme cases. However, at those energies where the curves indicate 
the maximum yield, the measurement of the curvatures was accurate to within 
4°. 'The discrepancy between experiment and theory noted in this region for 
lighter elements therefore seems real. 

In one aspect the present experiment confirms that of Alichanow, Alichanian 
and Kosodaew (1936), who observed the pair-formation in Pb by 2-2 Mev. y-rays 
from RaC and found that the maximum of the energy distribution curve of the 
positrons was at about 1-38 Mev. 


(ii) The Difference of Energy of Positron and Electron 


Approximate theoretical values for E,, —E_can be calculated from the formula 
given by Bethe and Heitler: 


paola eg 


where Z is the atomic number. 
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Table 2 sets out the experimental results together with those given by the 
theory. 


Table 2. E,—E_ 
Element Pb Pt Sn Ag Cu 


Theory 0-6 0:57 0:36 0-34 0-21 
Experiment 0-41+0:06 0:38+0-06 0-21+0-05 0-19-+0-05 0:10+0-04 


The experimental results appear to be consistently lower than the theoreticab 
values, but whether this discrepancy is real or whether it arises from approximation 
in the theory for the energy under consideration is not clear. The above formula 
is precise in requiring that the difference of the average energy should be propor- 
tional to the atomic number. 


A graph has been plotted to show the variation of FE, oie against Z. 


Figure 2. 


Figure 2 shows the graphs of both theoretical and experimental results. 
The theoretical results show a straight line with a slope h=0-0073. The experi- 
mental curve is almost a straight line having a slope h =0-0063 + 0-001. 


(iii) The Angular Distribution 


According to Bethe and Heitler, the fraction of the particles projected between. 
the angles @ and @ + dé to the incident photon is given by 


6d0 
(0? +02) 


where © =m,c?/E, where E is the total energy of either of the projected particles, 
i.e. the sum of the kinetic and rest energies. The maximum of the theoretical 
distribution of either positrons or electrons for y-rays of 2-2 Mev. when plotted 
against the direction of the incident photon lies at about 15°. From the experi- 
mental data it is found that the angular distribution of both positrons and electrons 
is the same in each element, and that the maximum of the curves of Pb, Pt, Sn, 
Ag and Cu lies at about 18°. Therefore only one curve has been used to illustrate 
the results. 
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Figure 3 shows the angular distribution for the positrons in Pb together with 
theoretical curve, given by the dotted line. The experimental curve follows 
the general shape of the theoretical distribution, and to this extent indicates fair 
agreement between the theoretical and experimental results. 
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Figure 3. 


The measurements of the angles between the directions of emission ot the: 
positrons and electrons are given in Table 3. 


® 


| Table 3 

Element Pb Pt Sn Ag Cu 
Average angle (°) between positron andelectron 31 30 30 28 28 
Average angle (°) between positron and y-ray 13 13 14 13 13 

Average angle (°) between electron and y-ray 18 V7 16 15 5 


The Table shows that the average angle between the positron and electron in 
the different elements is approximately the same, and is about 30°. ‘The measure- 
ments of the angles are accurate to within 3°. The mean of the average angles 
between the electron and y-ray is 16°, and that of the positron and y-ray is 13°. 

Klarmann and Bothe (1936) investigated the problem of pair-production in 
krypton and xenon by ThC” y-rays of energy 2:6 Mev. and found that the average 
angle between the positron and electron increased with the atomic number. 
Thus for krypton it was 19° and for xenon 42°, but the number of pairs studied 
was only 13. The present observation does not confirm this, but agrees with 
Simons and Zuber (1937), who found that the average angle in argon was 30° 
and in iodine 32°. The photograph (see Plate) shows the formation of a pair 
in copper by a y-ray of energy 2:2 Mev. 


§4. CONCLUSION 


Lawson (1949) has discussed the limitations of the theory of Bethe and Heitler, 
and it appears that this theory is not suitable for precise comparison with the 
present experimental results. For y-rays of energy less than 3Mev. Lawson 
states that the calculation of Jaeger and Hulme gives a correction of opposite 
sign to that required when comparing the experimental results with the theory 
of Bethe and Heitler. Since with heavy elements, however, the lack of agreement 
between theory and experiment possibly arises from the inapplicability of the 
Born approximation, while for light elements it may arise from the use of the 
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Fermi-Thomas atom, satisfactory comparison between experiment and any 
present theory is hardly possible. Bethe and Heitler’s theory predicts that the 
angular distribution should be the same for both electrons and positrons. ‘The 
experimental results show that the angle between the direction of emission 
of positron and incident y-ray is smaller than the angle between the electron 
and y-ray, but the difference observed falls within the error involved in the 
measurements and hence the agreement between the theory and experiment 
in this connection may be satisfactory. However, the most probable angle of 
ejection of either positron or electron for all elements is somewhat higher, as shown 
in Figure 3, than the theory indicates. In addition, the difference of average 
energy between the positron and electron is consistently lower than that to be 
expected from the theory. The change in the position of the maximum of the 
energy spectrum of the positrons with the element studied is quite distinct, 
and cannot be explained within the limits of probable error. 
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Structure and Diamagnetic Anisotropy of 
P-Benzoquinodimethane in Connection with those of 
P-Benzoquinone 


By M. G. EVANS; J. pr HEER{ anp J. GERGELY f, 
Communicated by C. A. Coulson ; MS. received 1st February 1949 


-ABSTRACT. ‘The interest in the diamagnetic properties of quinodimethane lies in the 
fact that this molecule can be considered as a prototype for quinone. The diamagnetic 
anisotropy, 4y, of p-benzoquinodimethane has been calculated by the molecular orbital 
“method following the treatment given by F. London in 1937. We find: 


A Xp-benzoquinodimethane ~ 0-34 Xbenzene 


‘Calculations of the diamagnetic anisotropy of quinone itself fail because the perturbation 
used here breaks down in this case, in which all the centres in the molecule are not identical. 

The result obtained for the quinodimethane can be discussed in connection with the 
observations by K. Lonsdale and K. S. Krishnan on quinone. We conclude that although 
-4y of quinone is much reduced compared with that for benzene, we need not postulate, 
to account for this observation, a complete or nearly complete localization of the zr-electrons 
in this structure. 


HOUGH Szwarc (1947) has only recently shown the presence of 
p-quinodimethane, 


CH,=< pci, 


in the products of the pyrolysis of p-xylene, this compound has been, for several 
reasons, studied theoretically since 1945 (Diatkina, Namiot, Syrkin 1945, Coulson, 
Craig, Maccoll, Pullman 1947, Pullman, Berthier, Pullman 1948). In the 
diagram below we show the free ealencies and bond-orders§ as calculated by 
Pullman e¢ al. (1948) by means of the molecular orbital method : 


-410 
=. 
758 yA ie 
TaN oe \ 
018 \\ 818 J 922 


For reasons which will be discussed below, we were interested in the dia- 
magnetic anisotropy, Ay, of this compound. A theoretical approach is the only 
‘possible one at present. Unfortunately so far no well-established general theory 
of diamagnetic anisotropy has been put forward. F. London, in 1937, developed 
a method for calculating Ay in certain aromatic systems with reference to Ay 
‘for benzene as standard. This method is essentially equivalent to a molecular 
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orbital treatment of the mobile z-electrons in these systems with the usual L.C.A.O. 
approximation.t— By choosing variation functions of a suitable type, the energies 
e! of the molecular orbitals, in a magnetic field of strength H in the direction 
perpendicular to the plane of the molecule, are obtained as the solutions of a 
secular equation of the general form 


| A 1 — €1 8 4| = 0. ecco ee (1) 


H }, reduces to the ordinary Coulomb integral, «, fork=1. Ayy= exp (27tfjq) 
= By if k and l are adjacent atoms (f is the ordinary resonance integral) and Hy, 
is zero otherwise. 5,,is the Kronecker symbol and the quantity f;,, above expresses 
in suitable units the magnetic flux through the triangle formed by the centre of 
the molecule and the bond between the two atoms k and J. In the case of the 
p-quinodimethane molecule an eighth-order equation is obtained, which can 
be reduced by means of group-theory (the only symmetry operation allowed in 
the presence of the magnetic field being an inversion with respect to the centre 
of the molecule) to the following two: 


Bo tl creed Wy AN) ae ay oy ome er | 
x E He x F a og 
P oa eee piped he pea Brice (26) 
0 nS * - % O nf «% 1% 
O nn ng * ) =, —Yi-* 


where x is written for («—e1)/8, the n*s are the complex conjugates of the ys, 
and the suffixes s and d indicate a ‘“‘ single” and a “‘ double”’ bond in thering of 
the molecule. Write 


mena =exp (if), sewn we (3) 
where f[(e/hc)H] is put equal to the area of the benzene ring. In this approxi- 
mate treatment we can neglect the small difference between the areas of the benzene 
and the “ quinoid”’ ring in quinodimethane. As f is a small quantity, we obtain 
from (2a), (2d): 

=Agt 2uile(arh toe? Gone ei cae oe (4a) 
Ax 2x l= = (ape. > eee (4) 


These equations, the solutions of which we will call «1, reduce to the ones we 
would have obtained from the ordinary secular equations for p-benzoquino-- 
dimethane if f were zero (H is zero). The solutions of the latter, which we will 
call x9, are known (Pullman, Berthier, Pullman 1948) and as f is a small quantity 
we obtain from (4a), (4 6) 


(xt — x0) /(nf)? = ee aie ee Tquiaiheeor (5a) 
(x1 — x)/(af)? = Pat ry pie (55) 


From (5a), (55) the x's and so the c's are readily calculated, and by summing these- 
over all the occupied orbitals we obtain the total mobile energy Ein the presence: 
of the magnetic field in terms of (zf)*. Finally, following London (1937), the: 


+ A summary of the approximations involved in this method is given by Coulson and Dewar- 
(1947) : i 
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diamagnetic anisotropy, Ay, is assumed to be equal to the diamagnetic suscepti- 
bility due to the ‘‘ motion” of the mobile electrons in the plane of the molecule, 
and thus is given by 


10E 
AX= — Ba Cer ONO 


The result of our calculation (which is expressed most simply in terms of Ayjonsene)» 
is EE encoquticd tacts i O27 Ax genzene 

This value is surprisingly low if we compare it with the experimental data 
for p-benzoquinone (such data lacking for the quinodimethane itself). Lonsdale 
(1939, see also Lonsdale and Krishnan 1936) reported a value for Ax, penzoquinone 
of 40-6 compared with 61-9 for benzene. ‘This relatively low diamagnetic aniso- 
tropy has been put forward by one of us (Evans 1946) as an argument to approximate 
a quinone as a molecule in which the z-electrons are completely localized in the 
“double” bonds. On the other hand McMurry’s investigations of the spectra 
of C=O compounds (1941) indicate that in a quinone z-electrons must occupy 
molecular orbitalsextending all throughthe molecule just as in the quinodimethanes.. 
Recently the authors have confirmed the latter point of view (Evans, Gergely,. 
de Heer 1949). We have shown that there still is an appreciable delocalization 
of z-electrons in a quinone structure even when we introduce in the secular 
equations appropriate parameters to take into account the extra binding energy 
of the C=O bond as compared with the C=C bond and the extra electronegativity 
of the O-atom as compared with that of the C-atom. But now the problem remained 
to account theoretically for the low diamagnetic anisotropy of these structures.. 
We tried to calculate the Ay for p-benzoquinone by introducing similar parameters 
to those mentioned above in the secular equations (2). In this way, however, we 
found a “‘negative diamagnetic anisotropy” (or a ‘‘ paramagnetic contribution” 
of the z-electrons), which of course has no meaning at all and at present only seems 
to show the limitations of the London method. 

The fact that this method breaks down for quinones is not surprising, as here 
a second-order effect is calculated by means of a first-order perturbation calculus. 
The fairly good results nevertheless obtained by London might be correlated 
with the fact that he only considered molecules which were described as 
“‘ crossable hydrocarbons”’ by Coulson and Rushbrooke (1940) and as “‘alternant. 
hydrocarbons” by Coulson and Longuet-Higgins (1947). These authors showed 
that only for these hydrocarbons is the molecular orbital method in its simplest 
form self-consistent. So, while this method is still useful to calculate quantities 
such as mobile energy, bond-orders and free valencies in both non-crossable 
hydrocarbons and hetero molecules (in the absence of an external field), it might 
easily give completely wrong results if used to calculate a second-order effect 
such as the diamagnetic anisotropy. 

As a result of this difficulty we had to confine ourselves to the calculation of 
Ay for p-benzoquinodimethane and to consider this molecule as the prototype 
of a quinone (as has been done before (Pullman, Berthier and Pullman 1948, 
Diatkina and Syrkin 1946)). As we saw, this computation gave us a very low 
Ax, though the relatively high resonance energy of 1-924 8 indicates an appreciable 
delocalization, and the bond-orders shown in the diagram above indicate that the 
ring in this molecule approaches more closely to a benzene ring than to a completely 
localized quinonoid structure. 
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Although, remembering the approximate character of the London method, 
-we do not attach any undue importance to the numerical value we have found, we 
tentatively suggest that a relatively small perturbation of the “pure benzene 
structure” in an aromatic compound can cause a rapid decrease in the diamagnetic 
-anisotropy. While a satisfactory explanation of this fact cannot be given at 
“present, we may conclude conversely that the low diamagnetic anisotropy found 
‘experimentally in p-benzoquinone need not be correlated at all with a complete, 
-or nearly complete, localization of the double bonds in this molecule. 


Note added in proof. It is easy to modify our calculations so as to include the 
-effect of an overlap integral S between the atomic wave functions on adjacent atoms. 
As is well known, the inclusion of this integral modifies the solutions, x9, of the 
unperturbed secular equation by a factor 1/(1—.Sx)). Brooks (1940) has shown 
that in a magnetic field the solutions («1 — x9) of equations (5a, b) above have to be 
“normalized” by 1/(1 —Sxp)?, the diamagnetic properties being far more sensitive 
to this non-orthogonality correction than the energy. Although with S=0-25 
the contributions of the individual electrons, <1, change appreciably, the total 
mobile energy £1 remains practically unchanged. 
Our final result now becomes 


SX penzéacihodine mane =0-25 SY henvene 
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The Computation of Wave Functions in Momentum Space—- 
I: The Helium Atom 


By R. MCWEENY * anp-C. A. COULSON t 
MS. received 19th April 1949 


ABSTRACT. The possibilities of obtaining precise momentum wave functions for atoms . 
by direct solution of the wave equation in momentum space are examined in some detail. 
An iterative method of approximating to the wave function is employed in computing the 
momentum distribution function for the helium atom. Although in this instance con- 
siderable accuracy is achieved, formidable difficulties arise in extending the calculations to 
more complicated atoms. 


§1. INTRODUCTION 


OR some time the importance of momentum wave functions has been 
appreciated: in particular, such functions are fundamental to any 
discussion of the Compton scattering of x rays (DuMond 1933) and the 

inelastic scattering of fast electrons (Hughes and Mann 1938). 

Hitherto the possibility of computing momentum wave functions directly, 
by solving the eigenvalue equation in the momentum representation, has usually 
been dismissed as impracticable and the required functions have been obtained 
by an application of transformation theory to the more familiar coordinate space: 
wave functions (Duncanson and Coulson 1945). Although this procedure is. 
satisfactory when the approximate solution is built up solely from hydrogen-like 
orbitals with appropriate screening constants, it proves impossible to transform 
to momentum space the most accurate wave functions in which electron cor-. 
relation is adequately represented, the integrals involved being quite insuperable. . 
This difficulty has prompted a further examination of the possibilities of a direct 
solution in momentum space. ‘The methods of calculation are essentially those - 
employed by Svartholm (1945) in the momentum space treatment of nuclear 
systems : a distinction arises, however, in that emphasis will be placed upon the: 
computation of the wave functions themselves rather than the associated energies ; 
for this purpose the iterative method is particularly suitable. In the present 
paper, the n-electron problem is formulated in momentum space and in 
the case of the helium atom the iterative method of solution is applied with. 
considerable success. In more general problems the “screening constant” ’ 
wave functions of Slater and others may be improved to some extent by iteration 
but there appears to be little prospect, at present, of obtaining from such improved | 
functions the necessary momentum distribution function. 


§2. FORMULATION OF THE PROBLEM 


For generality we consider a system of N particles, the first N—1 being electrons 
of mass m and charge —e, while the Nth is a nucleus of mass M and charge + Ze. 
Since it is impossible to avoid a certain amount of transformation theory a notation 
for matrix elements and transformation functions must be introduced : the Dirac 
notation seems most convenient for the formal derivation of the eigenvalue 


* University College, Oxford, now at King’s College, University of Durham. 
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equation. For practical purposes however, a familiarity with representation 
theory is not essential since we may begin our calculations starting with equation 
(8). In what follows we shall use either of the following equivalent notations for 
the momentum wave function and its complex conjugate : 


«Pix Pry + Rie -Pnx|)=P(Pix- gifts - Pz) 
(| Dis Pig Pu = 2 “(Piz Piva) 


The set of commuting observables p,,,....pPy, may suggest that @ is essentially 
a function of the Cartesian components of momentum : the choice of coordinates 
in which the wave function is expressed is, of course, arbitrary and the dependence 
of ® on the momenta may be abbreviated to ®=Q(p,.... py) i.e. a function of 
N momenta, each momentum being dependent upon three variables. In many 
cases the spherical polar coordinates (p,0,¢) of each momentum prove most 
convenient. Often the whole set of momenta will be abbreviated to p. The 
formalism should be clear from the following equivalences : 


{p|)=0(p)=O(piPpe2---- Py)=P(hr 91, $1 Pa» 9a $2, ---- by) 
dp=dp,dp,....dpy=p7sin 0, dp, d0, dd, p3....ddby. 
In momentum space the eigenvalue equation for the energy becomes 


i 5 1 N-1 5 
an Px t a ig pj |<p|)+ 
‘The representative (p|V|p’) of the potential energy operator is most expediently 
found by transformation theory. In coordinate space 


jal N-1 
Vi 2 Vin(|8in |) + » V5x( [Pix )s ee (2) 
Jal j<k 


jx, fjy being the vector separations of particles (j,k), (7, N). A single term may 
be transformed to momentum space according to 


{Pi-- +» Py [Vj Peps = eer | Cerin eee (Py ry )V(|r;—Ti|) 
Kb diy f,| Pine iy Pee eee ey 


‘where r,....fy are the position vectors of the N particles, and dr,....dry are 
volume elements in three dimensions. ‘The transformation functions are 


(p|r)=/-" exp{-ip.r/h}, <r| p) =k?" exp {ip.r/h}. 
The multiple integral (3) may be separated into a product of integrals over the 
3-dimensional subspaces and finally an introduction of relative coordinates gives 


(cf. Svartholm 1945) 
(Pi-- - Py |V. jh| Pie -- Py) 
=5(p;— Pj} + Pi — Pe Pi) 113(p Dice Di)ak Pilon) Be (4) 
p;) =h- Sexe ( =F pyar/h vs (le dese (5) 


With attractions of Coulombic form the appropriate integral is easily evaluated *, 
giving 


where (py |V 5x 


1 
*h | p;— p; |? 


* A convergence factor must be employed. 


Pili (| Ps) = — 
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Substitution of this result in (1) gives an integral equation for the eigenvalues 
and eigenfunctions. We confine ourselves to solutions in which there is no 
resultant translatory motion of the whole system by setting 


(Pi Pay -+++ P|) =8(Pi t+ Pe +--+. + Py) (Pi P2- ++» Pw-i|)- 
Putting p2= —2mE, R=m/M, A=m/n*h we find finally 


N=1 N-1 
| p3+2R Pia Pasa) aes #3 |p... Psd) 
ee = 
N—=1 
ar | = | <ps|Mawl Bi) ¢P)(Pi---- Pie Pru) 
N= 


N-1 ¢ 
+ = | <Pi| Viel Bi dP} Pr---- Phe. Pj+Pp-Pj---- py-al) | 
) y 4 


The accuracy of solution does not, in most cases, warrant the inclusion of R: 
iteration with a finite value of R would introduce a correction for the motion of 
the nucleus. We shall regard R as vanishingly small, and, with a slight change of 
variables, write the eigenvalue equation for an atomic system, for which the forces 
are Coulombic, in the following form, where x= N —1 is the number of electrons 
in the atom: 


BL n d 
Ep}(pi.---pal=al ZB | pe PLe P;—P---- Pnl> 
j=l1- 


n d SS 
= 2 [ie (pi--- y= P--- Pe FP» Pal) ee ee (8) 
j<z/ | P| 


§3. METHOD OF SOLUTION 
The integral equation (8) may be written quite generally as 


T(p)<p|)=A<p|V |p’ dap'{p'I). vee (9) 
In momentum space 7(p) is always algebraic, being in fact ahomogeneous quad- 
ratic form in the momenta, 

T(PiP2+-++ Pn) =(Po+Pi+P2+ ----Pr)- 
Solutions of (9) usually exist for a discrete set of negative values of EF (positive 
values of p%) : we are primarily interested in the (algebraically) least of these. 
For the purposes of solution it is convenient to regard A as the eigenvalue parameter : 
ultimately we determine an expression for an upper bound to the least eigenvalue 
° in terms of the constants of the problem and the unknown #3 : on giving X its 
true value 1/7? a corresponding value of p may be calculated. It is easily seen 
that the value of F obtained in this way is an upper bound to the least energy 
eigenvalue. 

The most powerful method of obtaining approximate solutions of an equation 
of this kind appears to be that contained in a theorem originally due to Kellogg 
(1922). Asimple proof is given by Svartholm (1945), who shows that the solution 
of (9) may be approached by iteration in the following way. 

Choosing an arbitrary initial function <(p|)° form the iterated functions 


: 1 
(pl) (pl)%.--+ according to the rule (p|)/'?= 7-5 | (piV’|p’>dp'¢p' |, 


and the integrals 


T;=J|p)T(p)dp<p|»”, 


W,=SJ<|p)dp<p|V |p’) dp'<p' |)’. 
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Next form the ratios 4;=7;/W,, Ajay=W;/Tji1G=0,1,....). Then the: 


sequence Ay>A,>A,>.... converges monotonically to the numerically least 
eigenvalue (A°) of (9). This result is dependent on the positive definiteness of the 
operators JT and V. 

For our purposes the trend of the iterated functions themselves is more impor- 
tant. It is not difficult to show that 


const. 


(Ply > Gag (19>, i>, 


where (p|0) is the eigenfunction corresponding to the eigenvalue A?. 
For, writing <p|)°==,(p|s)<s|)°, where |s) is the sth eigenstate of (9), 
and making use of the definition 7(p){p|s)=A°J(p|V |p’) dp’<p’|s), we have 


1 1 , , , (1) Sy 1 0 
(pl =2 775 | (PIV Ip) ap'<p'|s)<s|)°== 55 (P|s)<s| 


: j 
Repetition then yields (p|)’= (3) <p|s)<s|)>°, and since A° is by hypothesis 
§ 


the numerically least eigenvalue, the desired result follows immediately. 

The rapidity of convergence of the process thus depends essentially on the 
separation of the eigenvalues \°A1....A*.. .. (higher eigenvalues of A corresponding 
to excited states of a hypothetical system rather than of the actual one) and is 
difficult to assess precisely. It appears, however, that with a carefully chosen 
initial function one iteration is often sufficient to produce a rather accurate approx- 
imation to the eigenfunction. 

Even at this stage an advantage peculiar to the iterative method is evident. 
Although in many cases neither 7, nor W, is calculable, the formal expression 
for the eigenfunction with which they are associated, i.e. the eigenfunction 
corresponding to the third approximation, A,, to the eigenvalue, has already 
occurred in the calculation of Wy. ‘The superiority of this approach over a varia- 
tional method, in which the approximate eigenvalue must first be evaluated and 
secondly minimized in order to assess the best values of the parameters, is apparent. 


§4. APPLICATION TO THE HELIUM ATOM 
(i) Choice of Initial Function. 


With n=2 the eigenvalue equation becomes, using atomic units, 


(25 + P1 + Pz) O( Pi P2) =ALZ(y( pi Po) + Jo( pr P2)) —Li2(PiPs)], «+++. (10), 
where 


Up) Apa d es 
I= | pe O(Pi— P. P2), h= | eR» P2— Pp) 


Pee ia tally tant 5 
ha= | pers P,Pe+p), A=1/77. 


In constructing a first approximation to the wave function for an atom it is 
customary to regard the mutual interaction of the electrons as a negligible 
perturbation. Neglect of all integrals such as J,, then permits a separation of 
(10), or more generally (8), and indicates an approximate solution of the product 
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form analogous to the Hartree type solution in coordinate space. When spin 
functions are included and the required antisymmetry is introduced © is most 
generally a linear superposition of determinantal functions such as 


$1(pi)a(1),  ¢1(pi)B(1), Po(pi)a(1).-.. 
P= | 4i(ps)x(2), b1(Pe)B(2), ta(Pa)a(2)....f rere (12) 


eereee ec eee eee eee oer ees eee eee eeeesoee 


where ¢,(p), ®,(p),.... are the various possible wave functions for one electron. 
in the field of the nucleus. 


In the ground state problem at present under consideration 


b1(P) = $2( P) = $13 (P) 


and for a Hamiltonian not including spin variables, the antisymmetric spin 
factor in the wave function may, of course, be omitted. We may then write 


O°( p, Po) = $1s(P1)- Pis( Pe): 


As in coordinate space the single particle functions are of simple analytical form: in 
order to introduce some flexibility into the functions we admit a variable parameter 
and write 


1 
aN EP 
where, using atomic units, @ is numerically equal to an “effective” nuclear 
charge. ‘There is no need to normalize the function at this stage. 

Since this function is simply the Fourier transform of the screening constant 
wave function in coordinate space, we expect, and in fact obtain, an identical 
variational value of the energy from a calculation of Ay. The evaluation of the 
integrals 7, W, (or their equivalents), has been carried through independently 
- by Copson (1943): it is worth noting that the calculation is not appreciably more 
difficult than the corresponding calculation in coordinate space. 

In view of our previous observations as to the relative merits of variation and 
iteration, (13) appears to be the most suitable initial function, further refinement 
being most conveniently introduced by iteration. 


Op; ps) = 


(ii) The Iterated Function and Variation Treatment. 


Using the initial function (13) the first iterated function is given by 


(p53 + Pi + £3)O1( Py Po) =[Z(L1( Pi P2) +Lo( Pi Ps) —Li2(PiPs)], .--.-- (14) 
where, from (11), we find easily 
1 1 7 1 
T,( Pi P2) = a (peta) (pray > LoCpi Ps) = a Grey peta’ 
dp 1 
Lolpsps)=| nippy Fare Le ee (15) 


W,, may now be calculated: 


Wo = | [©°*(p5 +07 +p3)'dp, dpo. 
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By performing the integrations over p, and p, first the direct evaluation of /,, 
is circumvented and we obtain 


W, (ZIG A re (16) 


ial 


‘The calculation of Ty is straightforward 


T= os =; (2+ a). me: (17) 


1 2a?+p% Po 1 
Fant * Z—5/16 possesses a least value 2/2 (£516) 


SO Ay= when a?=93/2. 
Now, according to (11) we must put this minimum value equal to 1/7?._ This 
gives 


pi22(Z—5/16) AU, Par Sey Wg eee (18) 


In this way the most suitable value of a is obtained. In the actual computation 
of the iterated function we shall therefore take a=2 —5/16 =1-6875 a.v. 

Before proceeding with the computation a further remark upon the choice of 
numerical parameters is necessary, for (16) involves not only a but 3. 

Now in practice we cannot make an infinite number of iterations and so are 
never able to compute exactly the appropriate A in the relation A=A(p?). This 
implies that we can never exactly infer the inverse relation p3=p3(A) and conse= 
quently can never compute the true value of 3 in the case of interest where A = 1/7. 
Fortunately, however, the iterated wave function shows only a weak dependence 
upon p? and the use of a nearly correct p% instead of the completely correct one 
has no significant effect on the resultant wave function. In the present case we 
shall use the accurately known value pj = 5-8074 a.v. 

The significance of the various members of the iterated function becomes 
apparent on giving f@ its variational value, a?/2. With this somewhat inferior 
choice we find 


I,(PiP2) +/2( Pi P2) a 7 1 
(p° + Pi +P3) a (pi +a?)*(p5 +a??? 


which is simply a maps of the initial function. Accordingly we rewrite (14) 
in the form 


Oppo) =1'(PiPs)—Lio(PiPs)s wee eee (19) 
( 2 
Tol DBe)= Gao 


Here J’(p, pe) in itself represents a slightly improved wave function when fp? is 
given its true value: in this case 


pi +p3 +20? 


I(p.p,)= 
Pipe a (IPT PP +e peta? 
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which is no longer a simple product, f(p,).f(p,). In this way the probability 
distribution of p, is made to depend on the instantaneous value of |p,| : a similar 
small correction arises trom J{,, as may be verified by expanding the denominator 
of the integrand in terms of cos(p,, p), cos(ps, p): all terms beyond the first are 
dependent on the mutual inclination of the electron momenta as well as on their 
magnitudes. . 


(iu) The Momentum Distribution Function. 


Although the function ®1(p,p,) is of some intrinsic interest it is of little 


practical value: since the electrons are indistinguishable the only observable 


quantity is the probability of any electron having a given momentum p. ‘This 


> latter function is obtained by integrating the probability density over the momen- 


tum coordinates of all electrons but one. In the present case it is accordingly 


Pp,)=[©*(pips)®\(pipadpe wae (21) 


Moreover, since we invariably conduct our experiments on an array of 
randomly oriented atoms or molecules rather than on a single system, a further 
integration over the solid angle w, about p, is necessary. The quantity of 
greatest importance in the present case is thus 


P(p\)= | {©*(p.p.)® (p,p2) GD dae |) oe (22) 


P(p,) dp, then measures the probability of observing a momentum with magnitude 
between p, and p, +dpy. 

When simple product wave functions are employed the integrations are easily 
performed since the integrals separate into a product of ‘‘ one-electron”’ integrals 
{cf. Duncanson and Coulson 1945). In particular, we record the normalized 
‘probability distribution corresponding to the initial variation function: 


320 pe 
PO eae 


With a better wave function, however, such as ©'(p, p,), a separation of this kind 
is not possible and the integration over all p,-space may constitute a major 


| difficulty. 


Neglecting terms which are of the second order of small quantities 


Pi piy=s| (P2221 Teldpey y0. Gch ae (23) 


where J’” denotes the approximation to I’ which results on taking p?=a?/2; 


this choice of p3 in J’ makes a quite negligible change in the small correction term 


| 2I’T;, and at the same time greatly simplifies the computation. 


It is of interest to plot the corrected distribution retaining only the major 
term in the integral. We find easily 


2522 [seep 21(2/32, Pr) OMB) (24) 
@ | (tay! ~ “(Gita® ~ (+e? 


[r 2d Peo — 
where 
22 dz 


+ 0 
ei e)= | rare a 
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The various integrals have been tabulated as functions of p, and the modified — 
distribution curve plotted : to the scale adopted in the Figure the curve is, however, 


l | l | 

0 0:2 04 0-6 08 10 1:2 1-4 1-6 1:8 2-0 pay 
p (A.U.) 

1. Initial approximation P°(p). 2. Corrected function. P1(p). 


almost indistinguishable from the first approximation. 
The contribution of the second term in the integrand of (23) is 


2 | I'Tiodp, 


-={ lle 1 1 dpdp, 
(pi +a?) 


(po+@)? "(p> + Pi +3) [PP (Pi + Pp)? +4" }{(P2— PP +a yr” 


We may write 
me al . K'(p)dp : 
eT od. PERE WER Le SE ey oe 25); 
| whe a. (+a [p|?{(pi + py? +47}? si 
where 
1 l dp, 
K'(p)= = Faas Oia CLs 
(P GEp PD Ge) (pi py rae 
It is convenient to evaluate first 
dee so dp» Do aie 
K0)= GEER Ta rape aA 28 
from which K’(p) may be obtained by differentiating with respect to c and 
subsequently letting c—a. 
Taking polar coordinates about p 


fe oe Sdp, sind dé 
Kio) = P24P2 
(P= J. Jo), GERANGES EGET pP +a! Dpps cost 


The angle integrations are easily accomplished and the denominator of the 
resultant integral may be factorized to give 


K(p) 


oa) 3 dz 
J 0 (& +%c)(z—1c)(x +1b)(z —1b)(z — p +ia)(z — p —ia)(z +p +ia)(z +p —ia) | 
ive 
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By symmetry it is sufficient to evaluate the integral between the limits +00, this 
being accomplished by calculus of residues on completing the contour by an 
infinite semicircle in the upper half-plane. 

We proceed directly to the calculation of K’(p), writing down the residues, 
differentiating with respect to c, and letting c>a. After some reduction we 
obtain 


’ ae Apt + Bp? +C 
K"(p) = Cy Feestoceserrcrses Pa dap LDP pA =n |: eget (28) 


~where 


A=(a—b), B=(11at—16a%b + 2a2b? +304), 
C =(10a® — 32a°b + 30a%b? — 10a2b* +268) , 
P=(a2 +b), Q=(a2—B2). 


The'second integration, that indicated in (25), is rather more involved: it may 


tbe similarly reduced to a contour integral resembling (27) but with additional 


poles. The final result is 


4 
|rt.dp.— sapecapl ht RetR ves (29) 
where 
2 — 2B=Aa’) , 4(16Aat —4Ba? + C)(a? =p?) 
1— EM LM 
, (6Aat —4Ba? + C)(p} = 260%? + S7a4) 
4a? L7M 
p, _ AX+BY+CZ 
2 Dp tay’ 
Re a 
“2 (a+ Dy2ab(p? +a) —(a®— BH)’ 
and 


L=p;+10a*pi + 9at, 

M = pi — 6a°p} — 7a", 

X= (pp +@?)°(p} + 7a°pi +27a%pj + 210°), 
Y =(p} +.a°)?(p$ — api —Satp{ — 32°), 

Z = (p3— 1202p? — 42a*pt — 44a%p? — 158). 


"Remembering b= \/(p2 + p7) = \/(2a? +p{) to the present approximation, we 
may tabulate the above expressions as functions of p,. Since the result is angle 


_ independent the integration over the solid angle w, yields only a factor 47. The 


final momentum distribution including this further correction is illustrated in 
the Figure. The curve indicates a considerable decrease in the mean momentum, 
as would be expected from the virial theorem. For in momentum space 


il oo 
i oa P(p) p dp 


and a shift of the peak of P(p) towards regions of lower momentum implies a 
shift of weight from large to small values of p, with a consequent reduction in Hike 


’ Normalization is not easily effected analytically and for purposes of comparison 


tthe curve is plotted to the same peak height as with the first approximation. 
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(iv) Possibilities of Further Applications. 


We conclude with an enquiry into the possibilities of extending the analysis. 
just developed for helium to other atomic systems. Since the eigenvalue equation 
whose solution we are seeking refers essentially to a monatomic system, and 
association of any kind is known to have a profound effect on the momentum 
distribution, the next simplest atom to be dealt with should, logically, be neon, 
intervening atoms in the periodic table being incapable of sustained independent 
existence. 


It is essential that the initial function should possess the completely anti- 
symmetric form demanded by the Pauli Principle, for otherwise iteration would 
lead to a ‘“‘collapse”’ of the initial functions into the lowest available orbitals.. 
This naturally introduces enormous complications but there is still no formal 
difficulty in writing down an iterated function. Thus we find 


Re thes Steen tie Bae 
(PEPER soe Py (05 Pie ewe eee 


where S denotes the appropriate antisymmetric sum of a number of terms con- 
taining only one-electron integrals (every one of which is integrable with a 
screening constant initial function) and S,, denotes a similar sum of “correlation 
integrals’ (such as I,, of the last section). 

As already noted, the leading term in the expansion of each correlation integral 
could be evaluated and with the inclusion of such terms a rather accurate wave 
function might be obtained. 

In this way correlation effects would be represented to a considerable extent, 
and the corrected wave function would no longer be a simple product of one- 
electron functions. Unfortunately, however, aserious difficulty is simultaneously 
introduced. No function is of very great practical utility unless it permits the 
computation of P(p): the suggested function after iteration does not lend itself 
to such a computation, the task of integrating over all the momenta but one 
proving a formidable obstacle to further progress. For that reason we have not 
proceeded with the calculation of the iterated function. 


O'(p,. ees * P,) 
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ABSTRACT. A new approach is made to the problem of computing accurate momentum: 
distribution functions for electrons in molecules. The usual molecular orbital wave 
function (linear combination of atomic orbitals) has already been employed for 
this purpose but has proved inadequate in attempts to explain the width of the 
Compton profile. Using the method of a previous paper, the l.c.a.o. type of molecular 
orbital is formally corrected by iteration. In most cases it seems likely that the corrected 
orbital will be associated with a greater spread of the momentum distribution and will 
thus lead to a broader Compton profile : the broadening should become more marked with 
increasing delocalization of the electron and may apparently be very considerable. A com- 
plete calculation has confirmed these conclusions in the case of the hydrogen molecule ion : 
a preliminary treatment of the neutral molecule is given. 


SSS INEROD UC ETON 


N Part If of this work the distribution of momentum in a simple monatomic 

system was computed directly by approximate solution of the eigenvalue 

integral equation in momentum space; certain features of the iterative method 
employed were found particularly suited to this purpose. 

The object of the present investigation is to determine whether similar methods 
may be employed in computing the momentum distribution of electrons in mole- 
cules. A complete calculation is at present confined to the simplest of all diatomic 
molecules, the hydrogen molecule ion, but certain general effects of aggregation 
upon individual atomic orbitals are predicted. Results suggest that the use of 
more accurate wave functions might give very much better agreement than hitherto 
between the computed and observed profiles for Compton scattering from mole- 
cules and solid scatterers. 


§2. THE MOTION OF AN ELECTRON ABOUT 7” CENTRES OF FORCE 
(1) Formulation of the Problem 


For generality let us consider the motion of an electron about x centres of force 
located at points with position vectors R,, R,,...R,. Let the position vector of 
the electron be r.. We shall require the momentum space representative of the 
potential energy operator: this may be obtained using the methods and notation 
of I. If V, be the potential function for the sth centre 


V=X,V,(|¥ ==R; |). 
A single term is transformed to momentum space according to 
(p| V.(/r—R,|)| Pp’) =4-° fexp{—tp—p’.r/h}. V.(|r —R,]) dr. 


* Now at King’s College, University of Durham. 
+ Page 509 of this Journal: hereafter referred to as I. 
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With a change of variable we easily find, using atomic units as in I, ‘ 
{p| V.(/r—R,|)| Pp’) =exp{—7p—p’-R,}.(p|V(/|)|P), +++: (1) 
where ¢p| V,(|1|)| p’) is simply the momentum space transform of V, with the sth 


‘centre regarded as origin. 
The eigenvalue equation may thus be written 


(?3+p°)6(p) = —2m | exp{—ip=p’-R,}<p| Vip’ >ap'd(p’). -----2) 


If we suppose the electron to be in the vicinity of one particular nucleus and assume 
that all other terms in the potential function are therefore negligible compared 
with, say, V,, it is easily verified that 

$(p)=exp(—7p-R,).f(P) eee (3) 
is a solution of (2) provided ¢,(p) satisfies 


(p35 +p)d(p) = — 2m | (pl V.\ Pp’) dp’d.(P’). 
‘The function ¢,(p) is thus simply an eigenfunction for the motion of the electron 
about the sth nucleus alone, situated at the origin. We may therefore call 
exp (—zp.R,) a “displacement operator’. Putting 
D,(p)=exp(—ip-R,), Dy “(p)=exp(+7p-R,), 
the displacement of a physical system from 0 to R induces the following trans- 
formation of associated operators and wave function 


<p|V|p’> > Dp) p!| Vi p’)Ds(p’); 
¢(p) > D,(p). A(P). 


The eigenvalue equation (2) may thus be written 


(93 +P?)6(p)= -2m 3 | DPX PV. p')De(P)ap'$(p')..-.--(2) 


A superposition wave function precisely equivalent to the molecular orbital in 
coordinate space may be formed from terms such as (3). We shall call such a 
function ¢°(p) and take it as our initial function, namely 


°(P) =01D,(P)- $1(P) + o2D2(P). da(P) + --- 0, Dn(P)-$:(P)- +--+ (4) 
,(p)...%,(p) are the momentum space transforms of atomic orbitals about 
centres 1, ... m, and will in general be approximated by hydrogen-like wave 


functions. o,...0, are best assessed by the usual molecular orbital treatment. 


(11) The Iterative Correction 


‘Taking (4) as the starting point of our treatment the first iterated function 4'(p) 
may Bag ten be written down. From (2’) 


BAD) = Gee gy BEe | DARK PIV PYD,(P') dP DUP) H(P!) 
—2m 


me TP (Gea ph o> o,l"D.(p), 


where 
[t= | {p| V,|p’)exp{—7p’.R,—R,}dp'd(p’). —.... .. (6) 
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Even at this stage inferences as to the nature of the corrections may be made. 
-A typical term, JD,(p), of the summand in (5) may be said to represent the 
“perturbing influence of the sth nucleus on the momentum distribution of an electron 
“moving in the field of the ¢th nucleus alone. Even with neglect of all terms 
-except those for which s =¢ a correction results, for, to this approximation 


#4) = Geos 2 Dip) | (PIVp'>dp's('). (7) 

But, by definition of on atomic orbital 4p) 
[(e8).+2*1h(p) = —2m | (pl Vip’) ap’d(p'), eee (8) 
se ia oe Spe Dp yon ep). uslnaii a wt (9) 


= po+P" 


Now (9) is simply a linear superposition molecular orbital in which the individual 
-atomic orbitals have been modified by a factor which is an angle-independent 
function of p. Moreover, we can at once predict the effect of this correction on the 

momentum distribution function, for (p2), and pg may be regarded as, at least, 

approximately known quantities. Quite generally we may suppose pj>(p%), 

-since we are concerned solely with the ground state in which the electron occupies 

a bonding molecular orbital. Each atomic orbital in (9) is thus multiplied by a 

monotonically increasing function of p: consequently the resultant momentum 

distribution is a less rapidly decreasing function of p than the initial approximation 

‘would suggest. To a first approximation the half-width of the momentum 

distribution function is in this way always increased by the perturbing effect of 

neighbouring nuclei. Since, generally speaking, the difference between (p?), and 
_pe becomes more marked as the binding of the electron to the tth centre becomes 
less, we may say that the increase in half-width of the distribution function over 
that predicted on a simple I.c.a.o. treatment (Coulson 1941) will become greater as 
localization of the electron at a particular centre decreases. ‘This result is in 
agreement with experimental results on Compton scattering. Although the width 
-of the Compton profile is not simply proportional to the width of the momentum 
distribution there is always a marked similarity between the two curves and it may 
at least be shown that an increase in width of one implies an increase in width of the 
other. The computed profile agrees very well with experimental curves for 
monatomic scatterers but for molecular and solid scatterers the computed width, 
using the l.c.a.o. approximation to the wave function, falls far short of the 
observed width (Duncanson and Coulson 1941, p. 417). Assuming the approxi- 
mate additivity of one-electron distributions (Coulson 1941), these considerations 
suggest that refinement of the wave functions may lead to much closer agreement. 

‘This result is substantiated in a subsequent section. 

The modification of the initial wave function so far discussed may be described 
-as a spherically symmetrical distortion of each atomic orbital in the initial linear 
-combination. ‘The modifications imposed by those terms so far neglected are 

‘generally smaller and introduce an angle-dependent distortion, or polarization, of 
-each atomic orbital, together with a ‘‘mixing”’ of the orbitals. Reference to (5) 
‘shows that we maystill regard the molecular orbital asa linear combination of atomic 
-orbitals provided &,o,/* be interpreted asa distorted atomic orbital about the sth 
-centre. In addition to the dominant correction already discussed there will then 
‘be a contribution from every other centre, t, dependent on three factors: (i) the 
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form of the orbital about the ¢th centre, (ii) the relative importance of this orbital, — 
measured by o;, (iii) the magnitude and direction of the position vector of ¢ relative 
to Ss. 

In this way each atomic orbital takes up, on correction, the symmetry of the 
whole molecule. It is difficult to assess generally the effect of this correction on 
the momentum distribution function; calculation shows that in the case of the 
hydrogen molecule ion it is to oppose the correction first discussed. ‘The polari- 
zation effect will, however, become less marked with increasing nuclear separation 
and with increasing symmetry of the configuration, features characteristic of the 
solid state. 

§3. THE HYDROGEN MOLECULE ION 

For a homonuclear system such as Hj with nuclei 1 and 2 at —R and +R, 
(2) and (4) assume particularly simple forms. With a nuclear charge Z and 
attractions of Coulombic form we find, using atomic units 


: 2Z Aap ef 
(o8-+23)6p) = [SP ap oie), (10) 


oop) =2c0s (pRB). We a 6 Sie (11) 
where ¢,(p) is an atomic orbital for motion about either centre and symmetry has 
imposed the condition o, = oy. 

The momentum distribution corresponding to a wave function (11) has been 
obtained elsewhere (Coulson 1941). For purposes of comparison the normalized 
curve is plotted in Figure 1 as Pp). 

As we have already noted (I) a variational method of approximating to the 
wave function, using a linear function system such as 


pV P.,, (cos @) 
. m, 2 (pF cate ; 
is often precluded by the fact that it necessitates a computation of the energy 
eigenvalue and this may be vastly more difficult in momentum space than in 
coordinate space. For this reason we adopt the iterative approach. 

The first iterated function, corresponding to (5), is most conveniently written 


$O-pop leh ee (12) 
where I,(p) =2cos(p.R) _ apr Vip). oo a ee (13) 
and = T.,.(p) =2 Re [exp (zp. R)] | por exp (= 2p’: R).¢3(p')s Sas (14) 


I(p) and J,(p) then introduce, respectively, the angle-independent and the- 
angle-dependent corrections to the initial atomic orbitals. 
Since ¢,(p) satisfies, by definition, an ae 


[(P0)1 + 2° )1¢ -5 (2 ip-pP p o(P’), Par (15) 


where Z’ is an effective nuclear charge and a os associated value of p?, the: 
dominant term in (12) is 

[(Po)1 +P?) 
es Pps) 


AD Ot ee ht (16) 


: 
: 
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In the present case, guided by the screening constant calculation in coordinate 
space, we shall take Z’ =1-228. The corresponding exact solution of (15) for the: 
1s state is 


$(p) = 


(16) may therefore be written 


fees WIC Dg y= Zane I aera (en) 


2a" 1 
ss PR ie 7, 
a cos (p ) Brae TP) ( Po) 
Comparison with (11) shows that the principal modification introduced by 
iteration is the replacement of the simple 1s atomic orbital by the function 


ee 1 
$i(P) 5 ee ee eee (18) 


In coordinate space a more satisfactory molecular orbital might therefore be a 
superposition of atomic orbitals such as 


exp (= ar) — exp (=4r) (or 2 sinh (b—a)r 


tf Y 


.exp {—#(a+ iyr}) > 


this being the Fourier transform of (18), rather than ot the simple 1s orbitals,. 
exp (—a7). 

Since (18) is not appreciably more complicated than (17) we shall accommodate 
this result in our treatment at the outset by employing ¢/(p) rather than ¢,(p) as an 
initial atomic orbital. A further complete iteration will now be carried through 
with the inclusion of angle-dependent corrections. As in I we shall take the 
best available values of the parameters, namely p} =2:20523 a.u.; R=1-:000 av. 


(i) Spherically Symmetrical Correction 


We deal first with the angle-dependent part of the iterated function. With 
the new initial orbital, 4{(p), it is possible to evaluate I,(p) by standard methods 
of the calculus of residues. We find 


Ae 4 b)p 
1 
I,(p)=2cos(p. MLO: Pp tan 1 at. ae. dg) 
The numerical magnitudes involved allow us to replace the tan? by the first 
term in its expansion with a maximum error of the order 4%. Since it would be 
presumptuous to regard this as significant we take simply 


4772 


I,(p) = cos (p.R) eS 


Ultimately, as already noted, the momentum distribution function P(p) is 
required. From(12) 
2dw 
P(e) = | arom 
(P9 +2")? 
The term in J? arises, as we have seen, from the partially corrected initial 
function; the second term in the integrand introduces a further correction whose 


Peo eke oo ee (21) 
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order of magnitude justifies the neglect of the remaining, second order, term. We 
are here concerned with the angle-independent term: a simple integration gives 


Pp) = fe Paes P| a. sin dd db = pom? chi 5 ahaa |} e+ i]. 


This curve is compared with the screening constant approximation in Figure 1, 
where it is scaled to the same peak height. 


(11) Polarization Correction 


Before the correction due to polarization of the initial atomic orbitals can be 
‘computed it is necessary to obtain a series expansion of the integral [,,(p) which 
cannot be evaluated in closed form. Writing (14) as 


[,(p) =2 Re exp (7p. R). /(p), 
‘we may take polar coordinates about p as axis and, with R=(R, 0, 0), p’ =(p’, #’, ¢’), 
‘write 
ipl ie 2 { ” exp {—27p’R(sin 6 sin & cos ¢’ + cos 6 cos 0’)\p” dp’ sin 0’ d’ dg’ 
(p? +p’ —2pp’ cos 6")(p’”” +.a?)(p’? + B?) 


The ¢’ Peers on is easily accomplished (Sommerfeld’s Integral), giving 


ie 2x | re exp { —27p’R cos 6’ cos 6} . J 9(2p’R sin 6’ sin 6)p'"dp’ sin 0’ de’ 

oso (p? +p — 2pp' cos 6')(p™ + a®)(p +b”) 

aie (23) 

The remaining angle integration is awkward but may be peas term by term 
on making the expansion 

1 i 
(PIP*— pp cos) PPA ' 
which is obtained from a binomial expansion by expressing the powers of cos 6’ 


in terms of the Legendre functions. _ The ,’s are, of course, functions ot p and p’. 
Subsequently the following result, due to Gegenbauer (1874) may be employed 


P,(cosW)+ pq PkcosP ya. 5) cake (24) 


| exp (zz cos 0 cos #).J,_,(zsin # sin #). Ci(cos 0). sin’** 6 d0 
0 


: (z ) (isin? Cx(cos) .J,.,(2). 


Assuming the validity of (24) and remembering P,(cos 6’) = C3(cos 6’), an appli- 
cation of this result with v=} gives finally 


[A cos (p. R)—Bsin(p-.R)|p2dp’ 
ij Eo ae ee A a ee 
lee a: a) Cr lea) ck i ae 
where 

A = P)(cos 8) .J,(2p'R) — pp P(cos 0) . J 5/9(2p’R) + 4 P,(cos 8). J (2p'’R)—.... 
and 

B= — 1,P,(cos 6) .S s)o(2p'R) + j43P3(cos 8) . Szj9(2p'R) 

— H5P5(cos 8) J i1j(2p'R) + . 
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Expressing the y’s explicitly as functions of p and p’, and the Bessel functions in 
terms of elementary functions we might complete the integration expressed in (25). 
Since we require ultimately P(p) this would not be a very profitable course; 
instead we compute directly 


2p? 
Py(0) = Getom |, | LAPdy(e)sin 6.d8, 


with R as polar axis. The ¢ integration is trivial on account of axial symmetry- 
and we find finally 


si 1672p? 47? 
Pl) = Geremi| GabaETpy [Cl -DOV ------ (26) 
where 
C(p) = 25 Con Cras n even, 
D(p) =, Cae rns n odd, 
and 


+1 
Con =| COSa( PRX ioe, (¢).ds, ar evyen, 
sah 


=| cos(pRx)sin(pRx).P,(x)dx, nm odd, 
eet 
Crn =(—1)"? i RY PP’) -In+3(2p'R)p” dp’ 
i Ll \PR] (p? +a)(p? +B)(p? +p)’ 


n+1) ° 7 \E unl P')In+3(2p'R ped ( 
—(-0em | (oe) Gernot oe 


The functions C,,, are most easily built up from 


nm even,. 


1 rl 
In='| cos (2p Rx) 2x" dx, ~ n even, =| sin(2pke\.xvax; sntodd, 
e 0 


which may all be derived from J,, by the recurrence relations 

(2pR)Io, nui =("+1)Ig,—cos(2pR), 7 even, 

(2pR)lo, n41 =Sin (2pR) —(n + 16. n, n odd. 
The integrals Cp, involve the w’s: these are easily obtained explicitly by a power 
series expansion of (p? +p’? —2pp’ cos 0’)-linacos@’ where a=2pp'/(p? + p’”)( <1). 
The cosines may be written in terms of the Legendre functions by recursion. 
All the integrals are ultimately reducible to sums of integrals of the type 


6 p’’ sin (2p’R) dp’ 
ese || Se oe 
ee ae 
‘r cos (2p’R) dp’ 
I SBS A ele es 
(Cie)= |" pang 
Since a rather large number of these integrals must be computed as functions of p 
it is gratifying to find that for values of s>2 numerical differentiation of J(rS/2), 
I(rC/2) affords a convenient method of tabulation: higher integrals rapidly become 
smaller and this treatment gives adequate accuracy. The basic integrals J(rS/2), 
I(rC/2) may be evaluated by calculus of residues and have been tabulated for 
r=1,2,....5. A five-point rule proves adequate for subsequent differentiation. 


r odd, 


yY even. 
J 
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The results of the complete calculation appear in Figure1. Comparison 
-of the final curve with that obtained by Duncanson ( 1941) by transformation of the 


08 0 
p (A.v.) 


Figure 1. 1. Initial approximation Pp). 2. Partially corrected function Ps(p). 
3. Corrected function P(p). 


-coordinate space function due to James (1935) shows that the trend towards 
regions of higher momentum is maintained with increasing accuracy of the wave 
-function. 


§4. EXTENSION TO THE NORMAL MOLECULE 


In view of the preceding results it seems likely that even in the case of the 
normal molecule a preliminary correction of the initial function might be achieved 
by taking account of the angle-independent distortion of the atomic orbitals. 
Some progress along these lines has been made and will now be indicated very 
briefly. 

The eigenvalue equation for the normal molecule is found to be 


1 cos R 
(Po +Pi + P2)P(PiP2) = | 22 { | BPS dp; 0(p;P2) 


cosp>—p;,-R_, \ 
+] eet een IDL : 
ree wl p2 (pips) are p PP PsP. *P) | 


As in coordinate space, two types of initial function suggest themselves. ‘Taking 
‘the product of two molecular orbitals such as 


$°(p) =exp (ip. R). ¢,(p) + exp (—7p-R).41(p), 


and omitting terms corresponding to ionic states gives the ‘‘electron-pair” 
‘function 


(pipe) = 24;(Pi)¢i(P2)COSp;—Po-R.  —.......... (29) 


‘This function is known to be slightly superior to the complete product function 
and will be made the starting point of our treatment. The momentum distri- 
bution corresponding to (29) is easily calculated (Coulson 1941). With hydrogen- 
like 1s functions (17) we adopt Wang’s value of the effective nuclear charge 


* 
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Z' =a=1-166 but accept more accurate values of Rand pj, namely R=0-7000 a.v., 


- pj =3-7774 a.v., since we wish to compute the wave function for the known con- 


figuration of the system. 
The iterated function ©'(p,p.) may be written 


©(p,p,) = BEES [{341(p2) + 184(p1)} + (124s(P2) + 136,(P)} — Lal, 


where H}=2cospi—p,-R[ “PL dip) 
[pi—Pr| 


esa 
Ty =2Re [exp (ip; +P2-R)] | a ein) 


ET Sia Si 

interchange of suffixes yielding J$, [2, and 

exp (—2zp.R) 
Ip|? 


To compute P it is necessary to square the expression (30), integrate over all 
Pp 1 y to sq Pp g 
P.-space and finally over all directions in p,-space. Denoting the three terms in 


Typ =2Re[exp (iP: —P:-R)] | dp. $3(Pi =P) 44(Pa +P). 


Figure 2. 1. Initial approximation P%(p). 2. Partially corrected function Ps(p). 


(30) by E%, E? and £,,, and remembering E” and £,, are small compared with E° it 
is sufficient to compute 


° dp, pido, 

P(pi) = lomne, pay 

Js (pot+pi +P) 
‘The second and final terms in the integrand contain finer details of the corrections 
due to polarization and to electron correlation respectively: it 1s, however, 
noteworthy that the first and dominant term, in itself, 1s parametrically related to 
these corrections through p%. For the present we shall evaluate only the major 
contribution to the momentum distribution, retaining the first term of the inte- 
grand alone. In terms of the integrals 


[282 + 2b8k? —2EE |. 


+p’ sin 2p'Rdp’ prep dp’ 
rem acre (2) =| pata 
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where b? =p? + p?, a closed analytical expression for P(p,) is obtained: 
prpy _ SPR] sin 2pR) (H1S/42, pr) , 2108/32, pr) , 118/22, ps)) 
(=e lam CoB Grea + bray * Grae 
1(2/42, ps), 21(2/32, bs) , (2/22, eee |. 
(pi+a’)? (pi +a)? (pi +a") 
The integrals have been tabulated as functions of p, and the result of this first 
correction of the momentum distribution appears in Figure 2. The broadening 
of the curve is not so marked as in the cases of the ion: that is, however, a natural 
consequence of electron correlation since the mutual repulsion of the electrons 
tends to localize each of them in the vicinity of its own nucleus, making the mole- 


cular orbitals of the last section a less appropriate description of the behaviour of 


each electron. 
A more complete treatment of the problem would involve a rather compre- 
hensive tabulation of integrals of the form (14), and is perhaps not easily justified. 


§5. CONCLUSION 


The analysis in this paper makes it clear that the method of iteration is not very 
conveniently applied to any specific molecule: when it is remembered that our 
aim has been to construct really precise wave functions this is hardly surprising. 


Nevertheless iteration does appear to be of some value in indicating the type of 


correction which must be made in order to approach more closely the exact 


solution of a problem. Unpublished calculations show that the method may be 


applied with some success to problems having a more qualitative significance, in 
particular to the correlation of electron velocities by interaction. Only in this 
type of investigation does there appear to be any prospect of making any further 
useful progress. 
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Electrons Ejected by Polarized Radiation 


It is well known that, theoretically, the azimuthal angular distribution of photo- 
electrons ejected by plane polarized radiation follows the “‘ cosine-squared’’ law (Auger 
and Perrin 1927). The number of electrons ejected with components between y and 
x+dx is 

INCY OX ~COSUN GN, P= nai hitn y avany icdstene (1) 
where x is the azimuthal angle with respect to the electric vector E. (x is the component 


of angle measured in the plane perpendicular to the original photon.) |The electron 
intensity shown in polar form in Figure (a) is thus greatest in the direction of E. Kirchner’s 


(a) 


mean experimental results (Kirchner 1926, 1927), shown by the dotted curve, are in close 
agreement with theory. 

Comparatively little attention has been given to the azimuthal distribution of recoil 
electrons ejected by plane polarized radiation. It is, however, known that in this case the 
emission is greatest in the longitudinal plane containing the magnetic vector H (y=7/2) 
and not in that containing E. ‘The reason for this is that the scattered radiation intensity 
is greatest for y=7/2 or 37/2, and on Compton’s theory it is assumed that momentum is 
conserved in a collision between an incident photon and a free electron. The azimuthal 
angles of a recoil electron and scattered photon thus always differ by 7, whatever their 
longitudinal angles may be. If the further assumption is made, that to each scattered 
photon there corresponds one, and only one, recoil electron, the azimuthal distributions 
of recoil electrons and scattered photons will be identical, although of course the 
longitudinal distributions will be quite different. 

According to the Klein—Nishina theory (Heitler 1944) the number of photons scattered 
at a longitudinal angle @ into the element dQ, and with a given polarization, is 


Vo Vo Vg 


2 
N(0)xdQ ~ (24) (2 a Uo —2+-400s*X } OMe ten eee (2) 


/ where X is the angle between the electric vectors of the incident and scattered radiations. 
Following Heitler (1944), consider the scattered radiation as composed of two components, 
with polarizations E’ at right angles. For the first, E’ is perpendicular to E, so that cos K=O. 
For the second, E’ is in a plane parallel to that containing the direction of the primary photon 


PROC. PHYS. SOC. LXII, 8—-A 36 


530 Letters to the Editor 


and its electric vector E. In this case, cos? X=1—cos? x sin?@. Taking the mean of these 
components, the number of photons ejected in the direction (x, 8) becomes, from (2) 


Vo ° Vo Vo 2 : 
N(x, 9)dQ~{ 4) (= +——2cos*yxsin?6)dQ. —....... (3) 
Vo Vo Vy 
In (2) and (3), vg=vo/[1+a«(1—cos 4)], and w=hyvo/mc?, the energy of the radiation in 
terms of mc’. 
The azimuthal distribution of the recoil electrons, which is the same as that of the 
scattered photons, is now obtained by integrating (3) over 6 
0 


2 
n(x) dx ~ dx | (2) (2 + “2 —2 cos? x sin? 6) sar CHOCO, ek (4) 


Vo Vo Vy 
This can be put in the form 
UX): dX eo (LA cos? x) air) po pee ee ereine (5) 

The value of A in (5) can be determined analytically from (4) for the case of ~=0 
(¥g=Vo). Its value is then 2/3. This means that the ratio of the electron emission in the 
direction H to that in the direction E is 3 to 1. Graphical integration for other cases shows 
that, as « increases, the value of A diminishes. For w=1, A=0.52 ; for ~=10, A=0.16. 
At very high energies the electrons are emitted with azimuthal symmetry by polarized 
radiation. ‘The polar distribution is shown in Figure (b) for ~a=0.274 (140 kev.). ‘The 
dotted curve, which is in close agreement with the theoretical, represents Kirchner’s mean 
results for recoil electrons ejected in air and argon by radiation of this energy. 
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Theory of Photoconductivity of Layers of 
Semiconducting Substances 


In a previous paper (Schwarz 1948) it has been suggested that adsorption of oxygen ions 
during the evaporation of semiconducting substances and after the deposit has been formed 
is an essential condition for the production of photoconductive layers. Based on the experi- 
mental results, a qualitative theory of the photoconductivity of these substances is proposed. 

During the evaporation of the substance in an electric discharge, several types of adsorp- 
tion forces are probably active : (i) Van der Waals forces. (ii) Electrostatic forces between 
the positively charged evaporating particles and the negatively charged oxygen ions. Both 
kinds of ions are produced by electron bombardment in the’ electric discharge. (iii) Forces 
arising from the surface energy of the evaporating particles. (iv) Forces involving valency 
due to unsaturated valencies or to a small excess of the metallic constituent of the compound 
as a consequence of slight dissociation during the evaporation. 

Each particle when deposited is thus surrounded by an insulating layer of oxygen ions 
which form a coherent film round the semiconductor particles. | This intercrystalline layer 
fulfils three functions : (a) it prevents the coagulation of the particles; (b) by the formation 
of a double layer with the negative side away from the particle it constitutes a potential 
barrier for the exit of electrons from the particle into the neighbouring particles; (c) it is 
the source and carrier of electrons which by thermal excitation and/or radiation are liberated, 
accelerated by the external field and form a part of the current. 

The transport of electrons along a continuous layer of adsorbed oxygen ions may be 
visualized as hole conductivity ; a doubly charged oxygen ion loses one electron to a neigh- 
bouring singly charged oxygen ion, thus becoming singly charged, and in turn takes up another 
electron, 


ome 
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A cataphoretic migration of the oxygen ions in an applied electric field, especially at 
higher temperatures, is also possible. ; 

‘There are many observations on semiconductors which indicate the presence of a surface 
layer with special properties and there is strong evidence for a conductivity along the surface 
different from the bulk conductivity (Leo 1931, 1932, Dubar 1938). It is tentatively 
suggested that the hole conductivity as discussed above is at least partly responsible for the 
observed surface conductivity and that the surface levels introduced by Tamm (1932) and 
investigated by others are, in the case of photoconductive layers, occupied by oxygen ions, 
adsorbed on the internal surfaces of the semiconductor particles. 

The magnitude of hole conductivity will depend on the number of holes and on their 
activation energy, and these in turn depend on the thickness of the adsorbed layer, on the 
bond strength between the semiconductor particles and the oxygen ions and on the metal 
excess within the particle, as will be discussed later. 

The adsorption of gas on the particles during the evaporation, o1 in general before the 
bond between the particles can be established, is essential for the formation of a barrier layer. 
No heat treatment will produce a sensitive layer if the substance is perfectly free from oxygen 
and if no oxygen is admitted and adsorbed by the substance while the bond between neigh- 
bouring particles is interrupted. 

For the formation of an effective barrier layer and a photosensitive deposit it is further- 
more essential that a small metal excess within the particle is produced either during the 
evaporation or after formation of the deposit by a suitable heat treatment. If by unsuitable 
heat treatment the oxygen ions are removed by desorption, diffusion into the lattice or oxide 
formation, the photoconductivity is destroyed. 

The sensitive layer then possesses at the same time an excess of metal and non-metal. 
This small metal excess fulfils the following functions: (i) By producing a positive charge of 
the particle it facilitates the adsorption of negatively charged oxygen ions. (ii) It acts as a 
donor of electrons producing electrical conductivity of the particle. The activation energy 
will decrease with increasing metal excess and may even become negative. (ili) The work 
function may decrease with increasing metal excess. (iv) New absorption bands may be 
produced. Further experiments have to be made to clarify points (iii) and (iv). 

The final barrier layer has most likely a complex structure consisting of a chemi-sorbed 
layer near the surface of the particle and physically adsorbed layers on top of it. 

The bond strength between the oxygen ions and the particle decreases with increasing 
layer thickness so that the ions of the outermost layer approach the state of free oxygen ions 
for which the electron affinity for the second electron is negative. Although this state is not 
reached for adsorbed oxygen ions it indicates that the activation energy for the production 
of holes will decrease with increasing barrier layer thickness. With increasing metal excess 
and increasing barrier layer thickness the wavelength threshold for the same substance 
should therefore shift to longer wavelengths. 

Under the influence of radiation absorbed in the insulating layer positive holes are 
created which render this layer conductive and exert attraction on the electrons within the 
particle. At the same time radiation absorbed by the semiconductor particle causes an 
increase of the number and the energy of these electrons, thus increasing the pressure on the 
barrier layer. The potential barrier breaks down, and electrons can pass between the 
particles, constituting the photocurrent. It follows that any quantum yield greater than 
unity is possible, and that under illumination the ratio of excess to defect conductivity should 
increase. 

The electrical and optical properties of the photosensitive deposit will consequently 
depend on the distribution, size and properties of the semiconductor particles, on the pro- 
perties of the adsorbed layer and on the interaction between them. 

It follows from the foregoing that the photoconductive layer represents: a ‘state of non- 
equilibrium for the ionic and electronic constituents and may be considered as a disperse 
system. ‘The similarity in colour for layers of different substances and the kind of colour 
suggest that the size of the particles is fairly uniform for different substances and within the 
colloidal range. 

If adsorption is the essential factor in the production of these films one should expect to 
find these layers catalytically active. The high catalytic activity of many sulphides and 
oxides is well known, and the activating treatment of these catalysts closely resembles the 
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heat treatment of the photoconductive layers aiming at the appearance of a new phase— 
the small metal excess in the case of the photoconductive layer—without making it the pre- 
dominant component (Hiitting 1935). Rectifier plates of CuzO showed catalytic activity 
(Ostwald and Erbring 1931) and the efficiency of rectification increased with the activity, 
indicating that adsorption of oxygen ions may be essential for rectification of Cu,O plates. 
My thanks are due to Hilger and Watts Ltd. for permission to publish this letter. 


Hilger and Watts Ltd., E. SCHWARZ. 
98 St. Pancras Way, 
Camden Road, 
London, N.W. 1. 
22nd June 1949, 
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A Text Book of Heat, by LERoy D. Wetp. 1st Edition. Pp. x+436. (New 
York: The Macmillan Co., 1948). 25s. 


This is described as a book for upper classmen, and has been based on the material 
used for 35 years as lecture notes in a junior-senior course in heat. ‘This does not mean 
that it is old-fashioned, for modern work has been incorporated wherever it is needed. 
‘The properties of liquid helium receive due attention, the totally enclosed absorption type 
of refrigerator is described, and the International ‘Temperature Scale properly treated. 

The general outlook of the book is good. It aims evidently at presenting the student 
with principles rather than with practice, and many teachers may feel that it gives too little 
of experimental details; whether he is telling the student how to measure heat quantities 
or temperatures, or describing the evidence which led to the acceptance of the first law of 
thermodynamics, the author gives clear explanations of the principles of the experiment, 
but very little about the practical dodges adopted or required. 

The scope of the book includes thermodynamics and statistical mechanics, but the latter 
suffers because quantum theory is always treated as something outside, to be superimposed 
on the classical theory. Apart from this, it is an excellent treatment. 

‘The chapters are all accompanied with well selected sets of exercises, and with references 
for further reading. ‘These show a distinctly American bias, which is, indeed, evident 
throughout the book. The attribution of the theory of light darts to Epstein may be an 
example of this, or may be a misprint. 

As a whole, the book is very sound, and may give some very useful ideas and hints to 
those engaged in teaching this subject at the intermediate or pass degree level. ToElaan 


The Theory of Solutions of High Polymers, by A. R. Mituer. Pp. vii+118. 
(Oxford: University Press, 1948). 12s. 6d. 


This book summarizes, in a short space, modern ideas on the properties and behaviour 
of long molecules in solution. It can be said immediately that, so far as chemists are 
concerned, only those with an unusual flair for the mathematical approach will find the 
subject matter readily comprehensible. This is by no means the fault of the author’s 
presentation, which ts extremely good, but is due solely to the method of approach necessary 
in dealing with a difficult subject. Physicists and mathematicians cannot do otherwise 
than benefit greatly from having recent work brought together in this volume. 
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The author discusses first the classical theories of solution developed by van’t Hoff, 
van Laer, Gibbs and Lewis, and outlines the assumptions necessary for the derivation of 
such expressions as Raoult’s Law. He then explains the nature of an eigenfunction and 
of the grand partition function, and their application, first to general cases and then to more 
specific solutions. 

The second chapter is concerned with a description of long chain molecules and the 
connection between structure and configurational energy. Consideration is then given 
to the evaluation of the ‘‘ combinatory factor ’’, using Bethe’s method and modifications 
thereof. The method involves a statistical examination of sites occupied by molecules 
in a quasi-crystalline array. It is treated in detail for dimer and trimer molecules. ‘There 
follows a discussion of free energy of mixing of polymer molecules and the derivation of 
general formulae. 

The section having perhaps greatest interest, at any rate from a practical standpoint, 
is that which compares theory with experimental findings, e.g. vapour pressure, molar 
free energies of mixing and heats of dilution of such materials as polystyrene and natural 
rubber. Because of the nature of the initial assumptions, absolute agreement is not to be 
expected, but it is interesting that theory and practice are sufficiently related to justify 
confidence in the postulated ideas. Even for mixtures having non-zero energy of mixing, 
the relationships are sufficiently close, although the general impression given is that much 
planned experimental work remains to be done. 

The last chapter provides illuminating criticisms of the theory, and shows how the relative 
sizes of the polymer segments and solvent molecules, and the self-coiling of polymers 
may affect results. Comments are also made on the validity of the quasi-crystalline model 
adopted, and it is suggested that further advances may be made by devising a somewhat 
more satisfactory physical model. 

As stated above, the mathematical character of the book will make reading difficult 
for all except specialists in the field. To these, the matter summarized can only be of the 
greatest value; other workers interested in high polymers generally will be stimulated by 
the new concepts discussed by the author. ‘The style is logical throughout, the format 
is of the highest quality, and the absence of minor errors provides an excellent example 
of the care taken by the author in compiling this work. N. J. L. MEGSON. 


Radiology Physics, by JOHN K. Rospertson. Pp. xvii+323. Second Edition. 
(London: Macmillan, 1948.) 22s. 


During the last few years Professor Robertson’s Radiology Physics has been widely 
recognized and used as a valuable introductory course for medical and pre-medical students 
and for those studying the physical principles underlying Medical Radiology. 

The first edition, published in 1940, has been amplified very effectively by the addition 
of sections on T'ransmutation of matter and neutrons, Artificial radioactivity, mass and 
energy, and Uranium fission and atomic piles. In addition, the opportunity has been taken 
to meet previous criticism by the inclusion of a good deal of further information on ultra- 
violet and infra-red radiations and their applications in medicine. 

In accordance with the trend of events Professor Robertson has also included some 
discussion of the absolute energy systems of dosimetry, which seem destined to play a more 
important réle in the correlation of physical phenomena with biological implications. ‘The 
book does, therefore, provide an. admirable introduction to the more recent developments of 
pure physics and their applications to Medical Radiology, and will continue to be of great 
service to those commencing a study of the subject. W. V. M. 


X-Ray Optics, by A. J. C. Witson. Pp. vii+127. Methuen’s Monographs on 
Physical Subjects. 1st Edition. (London: Methuen, 1949.) 6s. 


The subject matter of this monograph is best indicated by its subtitle, “‘ The diffraction 
of x-rays by finite and imperfect crystals’. An opening chapter deals with crystal structure 
and the Laue diffraction theory and is followed by a discussion of the reciprocal lattice. 
This is used in the main part of the book, which is devoted to the structure of the diffraction 
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pattern from crystals with various types of imperfections, e.g. layer lattices with planes of a 
‘““ wrong ”’ type, alloys in a ‘“ disordered” state, distorted crystals. 

The mathematical treatment is kept as simple as possible. ‘The complex form for a 
wave is avoided throughout most of the book, and the reciprocal lattice is treated without 
vector algebra. It seems to the reviewer that the author is unduly pessimistic as to the 
mathematical equipment of those physics students who are likely to want to read his book. 
It is, however, a welcome sign that the need for the reciprocal lattice has at last been conceded 
by English writers on x-ray diffraction. It is to be hoped that the use of vectors will follow 
in the not too distant future. 

A somewhat disturbing feature in a lively and otherwise well written book is the curious 
terminology employed. No one will deny that physics is more important than phraseology, 
but is there really any need for the clumsy, misshapen and misleading expression “ three 
dimensional mistakes’, or, for that matter, for ‘“‘ frozen side bands ’’—another term 
favoured by the Cambridge School? M. BLACKMAN. 


CORRIGENDUM 


“The Establishment of the Absolute Scale of ‘Temperature below 1°k.”, by 
A. H. Cooke (Proc. Phys. Soc. A, 1949, 62, 269). 


Line 4, p. 274, should read 
Zo=2{1+2 exp (—6/RT)} instead of 
Z.=2{2--exp (—d/RT)}. 


Line 32, p. 276, should read 


4a le : 
X=Xo (: Be 3 Xo i 12ux3) instead of 


4a ; 2 
X=Xo (1 3 Xo! 12x) : 
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ABSTRACTS FOR SECTION B 


The Effect of the Surrounding Atmosphere on the Delayed Fracture of Glass, by 
C. Gurney and S. PEARSON. 


ABSTRACT. Round soda-lime-silica glass rods were subjected to a series of given bending 
moments, and times to fracture were recorded. Experiments were made in air and in vacuo 
of 10-1 and 10° mm. of mercury; and the effects of heating im vacuo to drive off absorbed 
atmospheric constituents were also investigated. For the same prior heat treatment the 
slopes of the curves of stress against time to fracture decrease with decrease of pressure. 
At 10-5 mm. Hg, after heat treatment, the curve was very flat. It is evident that the main 
cause of delayed fracture is attack of the glass by atmospheric constituents. Experiments 
in which glass was subjected to air from which water vapour and carbon dioxide were 
removed separately and together showed that both these substances cause delayed fracture. 
Other constituents of the atmosphere seem relatively ineffective. These conclusions are in 
agreement with the work of Preston and his collaborators. It is further concluded that these 
atmospheric constituents can be supplied by capillary liquid contained in surface cracks 
and by migration from surface layers as well as by direct attack from the gaseous phase. 


The Electrical Properties of Bismuth Oxide, by R. MANSFIELD. 


ABSTRACT. ‘The conductivity o and thermo-electric power dE/dT of bismuth oxide is 
measured over a temperature range of 680° to 150° c., and for variations of oxygen pressures 
from 76 to 10-4cm. Hg. The results above 340°c. and at atmospheric pressure are well 
represented by the formulae o= A, exp (—%/2T) and below 340° c. by c= A, exp (—«/T), in 
which the values of the constants A,, A, and « are very similar for two different specimens. 
The variation of o with oxygen pressure above 500° c. and for pressures down to 1 mm. Hg 
obeys the formula ao=k(Po,)*. ‘The thermo-electric power results may be summarized by 
the equation dE/dT=a-+b/T, where a and 6 are constants, the positive sign of dE/dT 
indicating that Bi,O; is a defect semiconductor. ‘The variations of dE/dT with oxygen 
pressure are small, less than 10%, except at pressures less than 10~* cm. Hg, when the sign 
of dE/dT changes at temperatures greater than 550°c. Under these conditions it is 
suggested that Bi,O3; becomes an excess or intrinsicsemiconductor. The results are discussed 
and good agreement is obtained between theory and experiment if itis assumed that normal 
Wagner and Schottky lattice defects become of importance at low temperatures when the 
concentration of free electrons is small, 
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The Mechanism of Dilatancy, by E. N. pa C. ANDRADE and J. W. Fox. 


ABSTRACT. The dilatancy in question is the change of overall volume produced by 
strain in an assembly of particles. ‘This has been demonstrated in the case of a mixture 
of sand and water by a photoelectric method of measuring the wetness of the surface. The 
general phenomenon has been experimentally studied in detail with a two-dimensional 
hexagonal array of uniform cylinders, part of the free surface of which is loaded with a rigid 
piston. An alternate increase and decrease of overall volume as the piston descends has been 
demonstrated and the mechanism of slip by which it comes about has been established. The 
elastic stress distribution in the array before slipping takes place has been calculated and 
shown to be a determining factor. The effect of friction between the component units has 
been studied. Experiments have also been carried out on three-dimensional irregular 
arrays of carbon shot and of sand. ; 


The Effect of Instantaneous Pre-Strain on the Character of Creep in Lead 
Polycrystals, by A. J. KENNEDY. 


ABSTRACT. ‘The extension-against-time curves of lead wires that have been subjected 
to rapid strain just before the experiment may be expressed by the Andrade creep equation 
1=1)(1-++ Bt*) exp (Kt), using the same constants as those which satisfy the creep of the metal 
under the same constant stress, but with the t value associated with 8 replaced by (t+ ty), 
where fy is a constant for a given experiment, its value increasing with increasing pre-strain. 
While for pre-strains less than 10°/, the value of « is unchanged by the pre-strain, for greater 
pre-strains the K-flow is also modified and appears to approach more nearly to a linear 
variation with stress. 


Observations of the Electrical Behaviour of Silicon Carbide Contacts, by E. W. J. 
MircuHey and R. W. SILLars. 


ABSTRACT. ‘The electrical behaviour of individual silicon carbide contacts has been 
examined using probes and high resistance electrometer circuits to measure the potential 
drop across a contact. Voltage-current curves were obtained and the well-known depar- 
tures from Ohm’s law at the contact were confirmed. Attempts were made to fit these 
experimental results and other data to the theoretical expressions for tunnelling through an 
insulating film, a hypothesis which does not appear to fit the facts very satisfactorily. 
Observations were made of the effect of heating the material to high temperatures in vacuo, 
of the properties of a newly formed surface and of the capacitance across the barrier layer. 
‘These indicate that the “ barrier ’’ has complicated properties which are not consistent with 
a simple tunnelling model. An alternative model is suggested. 


Temperature Measurements of Flames containing Incandescent Particles, by H. G. 
WoLruarp and W. G. PARKER. 


ABSTRACT. Colour temperatures measured on photoflash powders are higher than 
expected from theoretical calculations. Determinations of the true temperature of 
stationary flames of aluminium flakes suspended in air have been made (a) by a line 
reversal method and (6) by measuring the absolute light intensity. The experimental 
values obtained were very close to the theoretical flame temperature of about 3,000° c., 
that is, the boiling point of the oxide. Anexplanation is offered for the difference between 
the true temperature and colour temperature for flames containing incandescent particles, 
the reason being the small size of the radiant particles, which introduce optical anomalies. 
It is concluded therefore to be unprofitable to measure the colour temperature at all as this 
has no relation to the true temperature. Magnesium behaved very similarly, the true 
temperature being the melting point of MgO, near 2,800° c., whereas the colour temperature 
is above 3,900°K. The two-colour method of Hottel and Broughton for temperature 
measurements in hydrocarbon flames is examined in the light of this experience, 
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